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Abstract 



We review the phenomenon of color coherence in quantum chromodynam- 
ics (QCD), its implications for hard and soft processes with nuclei, and its 
experimental manifestations. The relation of factorization theorems in QCD 
with color coherence phenomena in deep inelastic scattering (DIS) and color 
coherence phenomena in hard exclusive processes is emphasized. Analyzing 
numerically the QCD evolution equation for conventional and skewed par- 
ton densities in nuclei, we study the onset of generalized color transparency 
and nuclear shadowing of the sea quark and gluon distributions in nuclei 
as well as related phenomena. Such novel results as the dependence of the 
effective coherence length on Q 2 and general trends of the QCD evolution 
are discussed. The limits of the applicability of the QCD evolution equation 
at small Bjorken x are estimated by comparing the inelastic quark-antiquark- 
and two gluon-nucleon (nucleus) cross sections, calculated within the DGLAP 
approximation, with the dynamical boundaries, which follow from the unitar- 
ity of the S matrix for purely QCD interactions. We also demonstrate that 
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principles of color coherence play an important role in the processes of soft 
diffraction off nuclei. 

I. INTRODUCTION 

The aim of this review article is to summarize fundamental principles of color coherent 
phenomena in quantum chromodynamics (QCD) and to apply these principles to the calcu- 
lation of cross sections of diffractive processes at high energy, where color coherence plays a 
crucial role. 

Our goal is threefold. Firstly, we make an attempt to give a comprehensive and peda- 
gogical overview of main ideas which constitute the concept of color coherence in QCD. In 
the preparation of this part, we have benefited from review articles |]],@], a talk ||, and a 
series of lectures f|. 

Secondly, we discuss hard exclusive diffractive processes. The recent proof of the fac- 
torization theorem for hard exclusive diffractive electroproduction of mesons || justifies the 
computation of the amplitude of this reaction in terms of the skewed parton densities 
In the present work, we analyze numerically the onset of generalized color transparency in 
hard exclusive diffractive processes and the phenomenon of nuclear shadowing in total cross 
sections in deep inelastic scattering (DIS) using the QCD evolution equation for the skewed 
and conventional nuclear parton densities. We also estimate the limits of the applicability 
of the DGLAP QCD evolution equation by comparing the calculated inelastic and elastic 
scattering cross sections of spatially small quark and gluon dipoles off the nucleon or nucleus 
with the dynamical boundaries, which follow from the unitarity of the S matrix for purely 
QCD interactions We demonstrate that studies of a new QCD regime, color opacity, 
may be feasible with future electron-nucleus colliders and in special hard processes in ep 
collisions at HERA. 

Thirdly, we review the soft diffraction of hadrons on light and heavy nuclei. We demon- 
strate that the principles of color coherence form the ground for the adequate phenomeno- 
logical description of the available data on hadron-nucleus coherent diffractive dissociation. 

The phenomenon of color coherence in QCD has two principal features: 1) an energetic 
projectile (hadron or photon) consists of coherent quark-gluon configurations of very different 
transverse spatial sizes; 2) a small-size quark-gluon configuration has reduced interactions 
with the target at not very small x. 

The coherent quark-gluon configurations, which constitute the energetic projectile, are 
called color fluctuations because QCD color occupies different volumes for the configurations 
of different transverse sizes. Phenomenon 2) is called color screening. Both aspects of 
color coherence have been proven in perturbative QCD (pQCD) and have been observed 
experimentally. Since configurations of various transverse sizes are present in the projectile 
at high energies, studies of color coherence phenomena provide an important relation between 
hard physics (small-size quark-gluon fluctuations), which is quite well understood in terms 
of QCD, and soft physics (large-size quark-gluon fluctuations), which is treated mostly 
phenomenologically. 

Property 1 is based on the observation that at high energies the longitudinal distance, 
essential for a given process, increases with the increase of energy. Originally, before QCD 
was known, such a behaviour of high-energy hadron collisions was suggested in P|-|n|] . The 
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longitudinal distance is called coherence length l c |T2|]. The word coherence reflects the 
attribute of any quantum field theory, which consists in the fact that particles are predom- 
inantly radiated at large and increasing with energy distances because of Lorentz slowing 
down of interactions. 

Property 1 is well understood in quantum electrodynamics ||13|| . Coherence length l c also 
exists in pQCD where it also increases with energy. It can be easily proven that this property 
of l c is a consequence of gauge invariance and the renormalizability of QCD. The convergence 
of integrals over "masses" of the intermediate states in the wave function of the incoming 
virtual photon follows from the absence of ultraviolet divergencies in the amplitudes, except 
where it is related to the renormalization of the coupling constant. However, in QCD this 
convergence is rather slow, which leads to the Q 2 dependence of the coherent phenomena, 
discussed in this review. This Q 2 dependence can also be understood as a result of the slow 
Q 2 dependence of the effective coherence length, introduced in the present paper. 

The large coherence length allows the presence of many quark-gluon coherent fluctuations 
of various transverse sizes in the energetic projectile. The latter results in a large probability 



of coherent diffractive processes, experimentally observed at HERA [ H| , and leading twist 



nuclear shadowing in DIS, experimentally observed for light nuclei JT^| and heavy ones [16| 



at CERN. Below we shall consider the relationship between the large l c and the experimental 
results in more detail. 

Let us consider diffractive DIS scattering of a photon with virtuality Q 2 , energy q , 
and high momentum q on a hadronic target T in the reaction 7* + T — > X + T 1 . The 
coherence length l c is introduced as the time (distance) for which the incoming virtual 
photon exists as quark-gluon fluctuations with the average mass squared (MJ-). It follows 
from the uncertainty principle that 

I — 1 _ 2q (1) 

C ^ Jq 2 + (M 2 ) - g <Mi) + Q 2 ' U 



where is the nucleon mass. In Eq. ([!]) one often uses the approximation (M|-) = Q 2 
and arrives at 

6 2m N x ' ^ ^ 

where x is the Bjorken variable |J. In Subsect. [11 B| we shall explain that since coherence of 
high energy processes is realized through diffractive processes, the parton model approxima- 
tion (M|) = Q 2 is no longer valid if one takes into account the effects of the QCD evolution. 
Then, more appropriate quantities are the effective coherence length l e J* defined by Eq. ([[]), 
where M\ is the typical mass in diffractive processes in DIS. 

Let the coherence length of Eq. ([!]) be larger than the typical size of the target 

lc > 2R T , (3) 

where Rt is the radius of the target. Then, not only are all quark-gluon fluctuations, 
satisfying Eq. (^), present in the photon wave function and interact coherently with the 
target T, but the target also has a large probability of remaining in the ground state. This 
fact leads to a large probability of coherent (when the nucleon or nuclear target is left intact 
during the interaction) diffractive processes on nucleons and nuclei. 
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In inclusive DIS on nuclei, when Eq. ([3D holds, the incoming virtual photon interacts 
coherently with several nucleons of the nucleus. This manifests itself in nuclear shadowing. 
The QCD-improved aligned jet model, which takes into account the interplay between soft 
and hard physics [fT^| , explicitly demonstrates that shadowing should be of the leading twist, 
i.e. only the logarithmic decrease of the effect of shadowing with increasing Q 2 is allowed. 
Note that at sufficiently small x the very concept of decomposition over twists might appear 
inapplicable to QCD. We refer the reader to Subsect. [IV for the relevant discussion. 

Property 2) of color coherence, color screening, has been proved in pQCD using gauge 
invariance, asymptotic freedom, and the QCD factorization theorem. The suppression of 
the interaction of small-size quark-gluon configurations of incoming hadrons and photons 
with hadronic targets leads to several distinctive experimental signatures of color coherence. 
In nuclear quasi-elastic reactions, where small-size configurations are produced as a result 
of high Q 2 two-body wide-angle scattering, the suppression of initial and final state inter- 
actions is called color transparency P"B| , |T5f . In QED, a similar phenomenon of "the electric 
neutrality" was theoretically predicted in as early as 1954 by H. Bethe and L. Maximon 
|20|j by considering coherent production of e + e~ pairs in photon-nucleus scattering within 
the eikonal approximation^. Note that one can justify the applicability of the eikonal ap- 
proximation in QED, while, in QCD, the validity of this approximation is questionable even 
in the perturbative regime because of "strong" gluon-gluon interactions. In hard diffractive 
processes, like the electroproduction of vector mesons off nucleons and nuclei, the interaction 
becomes strong at sufficiently small x, x < 0.01, and one reaches the regime of perturba- 
tive color opacity In this regime, the suppression in the amplitude is explained by 
nuclear shadowing of gluons and, hence, the factorization theorem is still valid. Therefore, 
the regime of perturbative color opacity can be also called the regime of generalized color 
transparency. We demonstrate that, in the limit of fixed x and Q 2 — > oo, the generalized 
color transparency is reduced to color transparency. At sufficiently small x, the interac- 
tion becomes strong, which leads to violation of the factorization theorem, and, possibly, to 
different nuclear shadowing in hard inclusive and hard exclusive processes. 

Perturbative QCD causes a particular energy and Q 2 behaviour of the differential cross 
section of the hard vector meson production. This behaviour has been verified experimentally 
(see discussion of relevant HERA experiments in JIIJ). Color transparency has recently been 
observed directly in the diffractive excitation of pions into two jets off nuclear targets |^3| 



in 



line with the predictions of |22| . A number of other processes such as J / ip photoproduction, 



leading hadron production in nuclear reactions using Drell-Yan lepton pairs as a trigger, 
and the production of leading nucleons in nucleus-nucleus collisions can be understood and 
confronted with the experimental data only if the coherent nature of the reactions is taken 
into account. 

This review is organized in the following way. Sect. |l| is concerned with cornerstones 
of the theoretical treatment of color coherent phenomena in QCD. These are the dynamics 
of small objects and the factorization theorems. We also discuss the relationship between 
diffraction on the proton in the reaction 7* + p — > X + p' and nuclear shadowing in DIS 
scattering 7* + A — > X. We critically examine the definition of the coherence length and 



1 We are indebted to S. Brodsky who drew our attention to |2C|. 
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present the Q 2 dependence of the effective coherence length l e J^ . We also review the aligned 
jet model of Bjorken and explain the improvements of this model which are to be done to 
make it conform with QCD . 



In Sect. pT| , the principles of color coherence in QCD are applied to exclusive hard 
diffractive processes. The main features of generalized color transparency are reviewed. We 
also discuss experiments, in which color coherence plays a crucial role. New reactions, where 
the two-gluon exchange is prohibited but color transparency still holds, are proposed. 

Sect. |Iy| quantitatively illustrates our conclusions concerning the transition from gener- 
alized color transparency to color transparency at fixed x and Q 2 — > oo, and concerning 
the behaviour of parton densities in nuclei by the virtue of analyzing the QCD evolution of 
skewed and conventional nuclear parton densities. A study of general trends of the QCD evo- 
lution in the x — Q-plane is presented. We also study the limits of applicability of the leading 
twist DGLAP evolution equation at small Bjorken x using the unitarity of the S'-matrix for 
purely QCD interactions as a guiding principle 0. 

In Sect. |V|, we give an overview of soft hadron diffraction on nuclei from the point of 
view of color coherence. We demonstrate that the formalism, which takes into account color 
fluctuations of hadronic projectiles explicitly, gives a universal description of the available 
data on diffractive dissociation and nuclear shadowing for light and heavy nuclei. 

II. COLOR COHERENCE IN DIFFRACTIVE PROCESSES 

The investigation of the role of color coherent phenomena in inclusive, exclusive and 
semi-exclusive high energy processes is one of the rapidly developing fields in QCD. Many 
calculations can be done in pQCD with the use of the dynamics of small-size quark-gluon 
configurations and the factorization theorems in a model-independent way. 

We begin this section by reviewing a deep relationship between nuclear shadowing in 
inclusive DIS on nuclei and DIS diffraction on protons, which is based on the factorization 



theorem for DIS diffraction 7* + p — > X + p' [23 . 



A. Leading twist nuclear shadowing and diffraction in DIS 

A deep connection between high energy diffraction and the phenomenon of nuclear shad- 



owing was first understood by V. Gribov [10,25]. On the qualitative level, nuclear shad- 



owing in DIS on the lightest nucleus, deuterium, arises from the interference between the 
amplitudes for diffractive scattering of the projectile off the proton and the neutron of the 
deuterium target. This is presented schematically in Fig. [IJ The application of the cutting 
rules of Abramovskii, Gribov, and Kancheli |26| explicitly demonstrates that the interfer- 



ence diagram decreases the total cross section for scattering, i.e. it leads to nuclear 
shadowing. 

One can use this result in order to calculate the modification of parton densities in 
nuclei at low values of Bjorken x. By comparing the QCD diagrams for hard diffraction 
and for nuclear shadowing due to scattering off two nucleons (see Fig. |l|), one can prove [27 



that, in the limit of low nuclear thickness, nuclear shadowing of nuclear parton densities 
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fj/A(x, Q 2 ) is unambiguously expressed through diffractive parton densities in the proton 

| p poo poo pXQ 

f j/A (x,Q 2 )/A = f j/N (x,Q 2 ) - - / d 2 b / dzi / dz 2 / dxjp 

2 J J —OD J 21 J X 

1 — rf 

x l+rj2 ff/N (A Q 2 , xp, t) | fc2=Q paO, Z!)p A (b, z 2 ) cos{xi P m N {z 1 - z 2 )) , (4) 

where fj/]\r(x,Q 2 ) is the conventional parton density in the proton; Pyi(r) is the nucleon 
density in the nucleus with atomic number A; xp is the fraction of the proton longitudinal 
momentum carried by the Pomeron; t is the momentum transfer, defined as —t = (k 2 + 
(xpmjv) 2 ) /(l — xjp), where k t is the transverse component of the momentum, transferred 
to the struck nucleon; (3 = x/xp] x ~ 0.02. In Eq. (Q), the factor (1 — r] 2 )/(l + r/ 2 ), 



where rj = —n/2d\n(^jf^ N )/d\n(l/x]p) = tt/2k accounts for the real part of the amplitude 
for the diffractive scattering P5[. The experimental value of the parameter k, defined as 



2k 



is rather small, k m 0.15. 



ff/N ((3,Q 2 ,X]p,t)^2 =0 OC X 

Considering deuteron, one can write a compact expression, which includes the t depen- 
dence of the electromagnetic form factor F2 H (g), where —t = q 2 

I pxq Y Tl 2 / \ 

f jf 2 H (x, Q 2 )/2 = f j/N (x, Q 2 ) - — J x dxjpdh^-^ ff /N (p, Q 2 , xjp, t) F 2H (4t) . (5) 

Similarly to Eq. (RI), one can easily generalize Eq. (H) in order to include the dependence of 
the diffractive amplitude on t ( see, e.g., P9fl ). 

The factorization theorem for DIS diffraction 7* + p —>■ X + p' states that the 
diffractive parton densities at fixed xjp and t are described by the leading twist DGLAP 
QCD evolution equation. Therefore, the nuclear shadowing, described by Eq. (f|), is also a 
leading twist effect and is governed by the DGLAP QCD evolution equation [[Z7j . 

Thus, Eq. ([D should be understood as the initial condition for the QCD evolution for the 
nuclear parton densities. Since Eq. (Q) is valid at low x, in particular at l e J$ > 2R A , where 
R A is the nuclear radius (see Eq. (^|)), the nuclear effects typical for larger x are absent in 
Eqs. (^) and (|Sp. A slight enhancement of gluon and valence quark parton densities in the 
nucleus at x ~ 0.1 should be added to Eqs. (|J) and (|5|) at some initial evolution scale Qq. 
At larger Q, this enhancement will contribute at much smaller x and will diminish the effect 
of nuclear shadowing. This simple argument shows that nuclear shadowing should decrease 
logarithmically as Q increases due to the QCD evolution due to two effects: decreasing l e Jf 
with increasing Q 2 and the growing contribution of the enhancement region x > 0.1. 

Eq. (|4]) can be used to study the amount of nuclear shadowing for quark and gluon parton 
densities separately. The analysis of []27| showed that: i) there is significant shadowing for 



quark parton densities with small theoretical uncertainties; ii) nuclear shadowing for gluon 
parton densities is larger by approximately a factor of 3 than for the quark parton densities^. 



2 Note that this result was derived in the approximation of the low nuclear thickness. Hence, it is 
valid for light nuclei only. For heavy nuclei, it is necessary to account for the interaction of 7* with 
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Fig. ^ presents the amount of nuclear shadowing at small x for the ratio of nuclear and 
nucleon gluon parton densities Ga/AGn and for the ratio of nuclear and nucleon quark 
parton densities qA/Aq^ for the nuclei of lead (Pb) and carbon (C). At the initial scale of 
the QCD evolution Qo=2 GeV, the nuclear shadowing corrections were calculated using Eq. 
(||) as well as taking into account the higher order rescattering terms. The dashed, dotted 
and solid curves correspond to Q=2, 5 and 10 GeV, respectively. One can see from Fig. |2| 
that at the initial QCD evolution scale Q=2 GeV, nuclear shadowing for gluons is larger by 
approximately a factor of 3 than for quarks. As explained above, due to the QCD evolution, 
nuclear shadowing for the ratio Ga/AGn becomes significantly reduced at Q=5 and 10 GeV 
and is compatible to nuclear shadowing for qA/Aq N . 

Note that, in many models of nuclear shadowing which do not use the information about 
gluon-induced diffraction, it is expected that gluons are shadowed less than quarks. Also, 
in a number of models, the interaction of small dipoles with several nucleons is considered 
in the eikonal approximation. However, by comparing Eq. (^) and the expression based on 
the eikonalization of the pQCD cross section for the interaction of a small-size dipole with a 
single nucleon (Eq. ([IB])), it is easy to check that the leading twist coupling of the small-size 
dipole to two nucleons, presented in Fig. [I], gives significantly larger shadowing than the 
higher twist eikonal diagrams, where four gluons are attached to a qq dipole. For example, 
at Q 2 =10 GeV 2 and x = 10~ 3 , nuclear shadowing, calculated using Eq. (|), is larger by a 
factor 3 than nuclear shadowing, calculated in the eikonal approximation. Hence, the eikonal 
approximation strongly underestimates the strength of nuclear shadowing for the interaction 
of small dipoles. In particular, for the total photoabsorption cross section of longitudinally 
polarized virtual photons <7l, Eq. (|J) corresponds to the leading twist shadowing for gluons, 
while calculations, using the eikonal approximation, predict higher twist nuclear shadowing, 
which is much smaller numerically, cf. |7j]. 



At the same time, it was recently shown |30| that a model, which incorporated the 
Gribov theory of nuclear shadowing and HERA diffractive data, successfully described the 
high precision NMC data on nuclear structure functions at low x [|T~5| , p~6[1 . 

In the laboratory reference frame, nuclear shadowing is controlled by Eq. (f|) and by the 
coherence length l c . In the next subsection we shall review the definition of l c and analyze 
its Q 2 dependence within the framework of the factorization theorem |24] and the QCD 
evolution equation. 



B. Q 2 dependence of the effective coherence length 

It is evident from Eq. (|]) that the leading contribution to nuclear shadowing in DIS 
on nuclear targets is proportional to the diffractive differential cross section at the minimal 
momentum transfer t m i n . Also, diffraction plays a dominant role in the generalized color 
transparency phenomenon. Thus, it is useful to introduce the notion of the effective co- 



3, 4, . . . nucleons of the target. Also, we have checked that the use of the diffractive parton densities 
obtained by Hautmann, Kunszt and Soper, based on the calculations in perturbative QCD at the 
initial scale Qo=1.5 GeV [84], lead to very similar numerical results. 
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herence length l e J^ . By definition l e J$ is the coherence length, which controls the onset of 
the nuclear shadowing and generalized color transparency phenomena. The definition of 
the coherence length l c , given by Eq. (^), is based on the parton model assumption that 
(M x ) ~ Q 2 . This approximation is valid only at moderate Q 2 and x ~ 10 -2 4- 10~ 3 . As Q 2 
increases, the ratio (M x )/Q 2 increases logarithmically at fixed x due to the QCD evolution. 
This effect becomes even more important at very small x due to the presence of a s \n(l/x) 
effects. 

According to the discussion in the Introduction, (M x ) should be understood as the 
average mass squared of the diffractively produced final state X in the single-diffractive 
scattering process 7* + T — > X + T'. Thus, in order to study the Q 2 dependence of l e J^ , 
one should determine the typical masses (M|-), which give the dominant contribution to the 
diffractive differential cross section da D /dt(t = t m i n ) of the process 7* +T — > X + T', where 

—tmin ~ (Xprn^) 2 . 

The diffractive differential cross section da D /dt(t = t min ) can be written as 
dv D ( _ \ d?a D _ 2ixa 2 , 2 \ r 02 P z>(4)/ n 2 

^ [t = tmin) = — 2 — [t = t min J = ~~^q^~ ^1 + (1 2/) J J dXjp P 2 \ X i Q 1 X Pi tmin) , 

(6) 



where F^\x, Q 2 , xp, t) is the diffractive structure function; y is the rapidity. For simplicity, 
the longitudinal structure function is neglected in Eq. (|6]). The Pomeron model of diffractive 
processes is supposed to be valid at small xp. We followed the prescription of Ref. |27]] and 
set the upper limit of integration in Eq. @ to 0.02. However, our numerical analysis is 
not too sensitive to the particular choice of the upper limit - a variation of the upper limit 
between 0.01 and 0.03 changes the result for l e J* by less than 25%. 



The factorization theorem for diffractive DIS scattering states that one can introduce 
the diffractive parton densities for fixed xjp and t, which satisfy the usual DGLAP evolution 
equation. The leading twist contribution to F 2 (x,Q 2 ,xjp,t) can be expressed through 
these parton densities. In practical applications, a further simplifying assumption is made - 
for sufficiently small xp < 0.02, where the vacuum pole exchange dominates, one can write 
F 2 D ^\x,Q 2 ,X]p } t) in the factorized form 

F 2 {i \x ) Q 2 ,x ]P) t) = fiP/vixjp^F^i^Q 2 ) , (7) 

where fjp/ p (xp,t) is the Pomeron flux; F 2 D ^ 2 \(3,Q 2 ) is the diffractive structure function; (3 
is the fraction of the Pomeron momentum carried by partons inside the Pomeron defined as 
(3 = x/xp. Substituting Eq. (|7|) in Eq. (^), one obtains 

da D . 2ira 2 / _ , 2 . Z" - 02 , „ , Sr ,D(2), x ^,2\ 

—^-{t = t min ) = — + -y) ) / dxp f ]P/p (x 1 p,t min )F 2 w ( — , Q ) 
dt xLj v Jx Xp 

= ^(i+(i- ! ,) 2 )/' ^/,^,w)if w w,<n- (8) 

xQ^ v Jx/om p z p 

In order to use Eq. (||) for studying l e J* and M\ as a function of Q 2 , let us notice that 
M\ is kinematically related to (3 as 



S 



This equation demonstrates that, at fixed by the kinematics of the reaction Q 2 and x, 
the averaging over diffractive masses Mx in Eq. ([]]) is equivalent to the averaging over j3. 
Substituting Eq. (||) into (|l|) we obtain 

l e J f = —(P)(Q 2 ,x) . (10) 
rriNX 

(P)(Q 2 ,x) is the typical value of (3, which dominates the integral in Eq. (|]). It depends on 
Q 2 and x. The Q 2 dependence of (j3)(Q 2 ,x) is governed by the DGLAP QCD evolution of 
the diffractive parton densities, which enter F^ 2 \(3,Q 2 ). The x dependence of (j3)(Q 2 ,x) 
comes from the x dependence of the lower limit of integration and the Pomeron flux in Eq. 

(!)• 

We propose the following ansatz for the evaluation of (j3)(Q 2 ,x). At each x and Q, 
(/3)(Q 2 , x) is given by the median of the integral in Eq. (g) 

of, l/p/pby^lAQ^/, -^fp/p{-3,tmin)F2 i2) W,Q 2 ) ■ (11) 

2 Jx/0.02 fJ z p Jx/0.02 p z p 

Thus, solving Eq. ( ITT]) for (j3)(Q 2 ,x) will give the desired Q 2 dependence of l e J* through 
Eq. ©• 

In our analysis of Eq. ([TT|), we used the leading order expression for the diffractive 
structure function F 2 D<y2 \f3,Q 2 ) 

F 2 D{2) (P,Q 2 )=Y,^PQa(P,Q 2 ) , (12) 

a 

where q a (P, Q 2 ) is the parton density of a quark with flavor a inside the Pomeron; e a is the 
charge of the quark. 

The Q 2 dependence of q a (/3,Q 2 ) was studied using the next-to-leading order (NLO) 
DGLAP QCD evolution equation. The recent parameterization of diffractive parton densi- 
ties of Ref. |31] was used as an input for the QCD evolution. We choose the fit D, which 



gives the best description of the diffractive DIS and photoproduction data from HERA. All 
fits of |H] correspond to the initial scale Qo=2 GeV. At this scale, the gluon density is 
non-zero, u(x) = d(x) = u(x) = d(x) ^ 0, and s(x) = s(x) = 0. As one performs the QCD 
evolution to larger Q > Q , the equality u(x) = d(x) = u(x) = d(x) remains and strange 
quarks s(x) = s(x) ^ are generated. 

The current HERA data does not cover the region of small f3 at small x, hence the fits 
of J3l| are not sensitive to the region where the triple Pomeron contribution dominates. 



In order to explore the sensitivity of the value of l z J* to this region, we have included the 



small-/? contribution, using the limiting fragmentation ansatz of Ref. [32]. In the limit of 
small j3 and at Q 2 = Qq, the multiplicity of leading hadrons in the target fragmentation 
region approximately coincides with that in photoproduction. This prediction is consistent 
with recent HERA data at xp > 0.05. 

In Fig. we present V^^m^x = (f3)(Q 2 ,x) as a function of Q at x=10~ 4 , 10~ 3 , and 
10~ 2 , calculated using Eqs. (0) and ( [TT]) with (dashed lines) and without (solid lines) 
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small-/? contribution. One can see that, at each x, the coherence length l e J$ is a decreasing 
function of Q. Also, one observes that the small-/? tail becomes increasingly important with 
decreasing x, especially at Q 2 ~ Qq. 

Using Fig. |3|, one can also study the x dependence of l e Jf at given Q 2 . While l e J^ increases 
as x decreases, this increase is slower than 1/x behaviour predicted by Eq. (0). Also, the 
difference between l e Jf for different x becomes larger as Q 2 increases. These effects are 
attributed to the QCD evolution of diffractive parton densities. 

Note that the definition of l e J^, which uses Eqs. fllOD and (|TTD, is less reliable at large 
Bjorken x, x > 0.01, since the diffractive cross sections are quite small at those x. 

We would like to stress that a decrease of l e J$ with increasing Q is a model-independent 
statement. It is a direct consequence of the factorization theorem in diffractive DIS and the 
QCD evolution of diffractive parton densities. As Q increases, the parton densities become 
softer, i.e. they shift towards lower (3 (or higher M x according to Eq. Hence, at fixed 
x, (/3)(Q 2 ,x) decreases as Q is increased, or, alternatively, (M|-) increases as Q is increased. 
Hence, by virtue of Eq. (P l e J> decreases with Q. 

Such a behaviour of the effective coherence length is supported by the experiment. The 
ZEUS data on single-diffraction |33| demonstrates that as Q 2 increases, diffractive states of 
large mass become increasingly important Hence, according to Eq. (|I|), the effective 

coherence length should decrease as Q 2 increases due to increasing (M|-) with Q 2 . 

In conclusion, we would like to point out that the present-day knowledge of diffractive 
quark and gluon parton densities enables one to evaluate the leading twist contribution to 
nuclear shadowing in the quark and gluon channels separately. In the nuclear rest frame, the 
quark channel corresponds to quark-antiquark fluctuations of the incoming virtual photon, 
while the gluon channel can be probed via a study of the quark-antiquark-gluon fluctu- 
ations of the incoming photon [J35| . The fit of Ref. |TT| suggests that the shapes of the 



/^-distributions of the diffractive gluon and quark densities at low Q 2 , Q 2 ~ 4 GeV 2 , are 
similar. This leads to the similar coherence lengths (up to the effects related to the 
region of small (3, which are somewhat different in the two cases). With an increase of Q 2 , 
larger scaling violations in the gluon channel lead to a faster decrease with Q 2 of l z J* for the 
gluon channel than for the quark channel. The detailed Q 2 dependence of l e J* for the gluon 
channel will be presented elsewhere. 



C. The dynamics of compact systems 

As discussed in the Introduction, at high energy, a projectile consists of coherent quark- 
gluons configurations of very different transverse spatial sizes. One can select rare processes 
where the small-size, or compact, quark-gluon fluctuations either dominate the process as 
in electroproduction of vector mesons by longitudinally polarized photons with high Q 2 
or are as important as the large-size configurations as in inclusive DIS or in diffractive 
processes initiated by transversely polarized photons. The cross section of the interaction 
of the small-size configuration with the hadronic target is calculable in pQCD since small 
distances correspond to large Q 2 , where pQCD is applicable. It was found that the small-size 
configurations have reduced interactions with the target. 

The cross section of the interaction of a compact qq fluctuation of the virtual photon with 
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a hadronic target T was derived in the leading a s log(Q 2 /A 2 ) approximation in [f?j,p^, |36| , j37l1 f| 



af (6 2 , x) = ^-b 2 \xG T (x, A/6 2 )] a s (\/b 2 ) , (13) 

where b = (b q — bg) is the transverse inter-quark distance within the wave function of the 
photon; a s (A/6 2 ) is the QCD coupling constant; Gt(x, X/b 2 ) is the gluon parton density of 
the target. This formula is a consequence of gauge invariance, asymptotic freedom, and the 
factorization theorem in QCD. The QCD evolution is properly taken into account through 
the gluon distribution of the target Gt(x, Q e ff = A/6 2 ). The scale factor A can be estimated 
from the analysis of the erj,(7*iV) cross section j7|. A slowly increases with the decrease of 
x due to the increase of the gluon density at small x. For x ~ 10~ 3 and Q 2 ~ 10 GeV 2 , 
A ~ 10. For a recent detailed study of an approach to inclusive DIS based on Eq. (|13|), see 

[sg. 

For a colorless dipole consisting of two gluons, the dipole cross section will be larger by 
a factor of 9/4, which is the ratio of the Casimir operators of the SU(3) C color group for 
the octet and triplet representations 0. Note that this observation is consistent with the 
analysis of DIS data in Ref. |?7J which demonstrated that the effective dipole cross section 
in the gluon channel is larger (by approximately a factor of 3) than in the quark channel. 

The gluon density of the target in Eq. (|T3|) is a feature that is absent in the two-gluon 
exchange models [|39|. The two-gluon description in [^j is incomplete because at small 



Bjorken x, a compact quark-gluon configuration interacts with the target by producing a 
soft self-interacting multi-gluon field. The effect of this soft gluon field is expressed by the 
gluon parton density of the target xGt(x, A/6 2 ). 

The proportionality of a^{b 2 ,x) to the gluon density xGt(x, A/6 2 ) has found many ex- 
perimental confirmations in inclusive DIS and in hard exclusive reactions on nucleon and 
nuclear targets. 

Let us consider the implications of Eq. (|T3|) for inclusive DIS scattering and for the 
famous approximate Bjorken scaling for the total photoabsorption cross section <Jtot{l* +P — > 
X)(Q 2 ,x). According to Eq. (|13|), the cross section for the interaction of a compact quark- 
gluon configurations with a target is proportional to the transverse size 6 2 oc 1/Q 2 of the 
compact configuration. It is important to note that this non-trivial result contradicts our 
experience from the models, which existed before QCD. For example, in the meson cloud 
models of the photon wave function, all meson configurations interact with the target with 
typical hadron-nucleon (nucleus) cross sections. In the case of DIS on nuclei at not extremely 
small x, such a behaviour leads to a to t oc lnx at x <C l/IR^m^ (see Eq. (||)), which 
contradicts the experimentally observed Bjorken scaling and causes a paradox. The paradox 
is known as the Gribov paradox. 

QCD explains the origin of the Gribov paradox. Indeed, in pQCD, the contribution of 
small-size configurations to 0-^(7* +p — > X)(Q 2 , x) behaves as 1/Q 2 according to Eq. (|13|). 
As a result, it has the same Q 2 dependence as the non-perturbative contribution of large- 
size configurations to o- to t{.l* +P — ► X)(Q 2 ,x), which is suppressed by a factor proportional 
to 1/Q 2 since the large-size configurations occupy a small fraction ~ k 2 /Q 2 of the phase 



3 This formula is also valid in the leading a s log x approximation. 
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volume. Therefore, the sum of these two contributions behaves approximately as 1/Q 2 at 
large Q 2 in accordance with the approximate Bjorken scaling. 

The QCD prediction of Eq. and its experimental confirmation both demonstrate 
that the existence of quark-gluon configurations, interacting weakly with hadron targets, is 
a fundamental property of QCD. 

Eq. (|13D substantiates the pre-QCD hypothesis of J. Bjorken, which became known as 
the aligned jet model [J40| . Within the framework of the parton model, Bjorken offered a 
solution of the Gribov paradox by suggesting that the interaction of the quark and anti- 
quark with large relative transverse momenta p t , which form the photon wave function, with 
hadronic targets is suppressed by some factor f(^ 2 /(p 2 + yU 2 )), where fi is a cut-off scale. 
Since large relative momenta correspond to small relative distances, the aligned jet model 
foresaw the suppression of the interaction of the spatially small quark- ant iquark dipole. So, 
in the aligned jet model, only large-size configurations take part in the interaction with the 
target. Since the available phase space of the aligned jet configurations is proportional to 
1/Q 2 , the scaling of (?tat{l* +P ~~ * ^){Q 2 i x ) is restored in the aligned jet model. 

In QCD, Bjorken's aligned jet model needs to be modified. Firstly, Bjorken scaling 
of (Jtot(l* + P — ► X){Q 2 ,x) is only approximate due to the QCD evolution with log(Q 2 ) 
of parton densities of the target. Secondly, as one can see from Eq. flT3|), the interaction 
of small configurations is suppressed but not negligible. Thirdly, as a result of the QCD 
evolution, the gluon parton density of the target becomes large at small x, which makes 
the cross section sizable even for relatively small b. Fourthly, one has to account for the 
suppression of on-mass-shell quark production and for the compensation of this suppression 
by the radiation of hard gluons ( see, e.g., |r?|J¥l| ). 

Leading twist nuclear shadowing predicted by the QCD-improved aligned jet model has 
been observed in inclusive DIS scattering on light and heavy nuclei, see ]IS],|TBJ and references 
therein. 



D. The factorization theorem for hard exclusive meson electroproduction off 

hadronic targets 

A crucial element of perturbative calculations in QCD are factorization theorems. Such 
theorems prove that amplitudes of certain hard processes can be computed as the convolution 
of the hard part, which can be computed in pQCD, and of the soft, non-perturbative, part, 
which contains information about parton distributions in the target. 

Recently J. Collins et al. || have proved the QCD factorization theorem for hard exclusive 
meson production in the reaction + T — > M + T", where M is a meson (p, n, K, r], . . .) 
and T and T' are baryon states of a fixed mass. This theorem is valid for hard processes 
initiated by a highly virtual longitudinally polarized photon on any hadron target at fixed x 
and Q 2 3> Aq C£) . The amplitudes of such processes are calculable as the convolution of the 
hard kernel Hij(Q 2 x\/x, Q 2 , z, /z 2 ), calculable in pQCD for a given minimal Fock component 
of the light-cone wave function of the meson wave function <j)j(z, /x 2 ), and the skewed parton 
densities of the target fi/ p {x\,X\ — x,t,fi 2 ) (see Fig. £| for the schematic representation of 
the amplitude for the process 72 + T — > M + T') 



12 



A = J2f dz J dx 1 f i/p (x 1 ,x 1 -x,t,^ 2 )H ij (Q 2 x 1 /x,Q 2 ,z,fi' 2 )(l)j 4 (z,[i 2 ) . (14) 

In Eq. (PI) , z is the fraction of the photon longitudinal momentum, carried by the quark 
or antiquark that constitutes the meson M; /i 2 is the factorization scale which is usually 
set to Q 2 ; x± and the longitudinal light-cone momentum fractions of the partons 

connecting the hard and soft parts of the amplitude A. Note that the skewed parton densities 
evolve as a function of x± and Q 2 for fixed x = Q 2 /2(p ■ q), while the conventional parton 
densities evolve as a function of x and Q 2 . In order to emphasize this difference, we keep a 
different notation for the skewedness parameter A = x\ — [x\ — x) = x. 

A similar QCD factorization theorem can be derived if M and T are baryons and T 1 is 
a meson. 

Applications of the factorization theorem give a practical possibility of calculating cross 
sections of hard exclusive processes in the same sense as total cross sections of inclusive DIS 
|5|-[7|]. The kernels for the QCD evolution of the skewed (nondiagonal, nonforward) parton 
densities have been calculated in the leading a s log(Q 2 /A 2 ) approximation in [P^ , pffil and 
in the next-to-leading a s log(Q 2 /A 2 ) approximation in |]44] , ^5| . Thus, the major difference 
between computing cross sections of hard diffractive processes and more familiar inclusive 
DIS processes consists in the necessity to use the skewed parton densities instead of the 
conventional parton densities. 

A natural assumption for the input skewed parton densities fi/ p (xi,xi — x,t,fi 2 ) at the 
initial scale /x = Qq of the QCD evolution and small x is that they should be taken equal 
to the usual parton densities at Q [(43]. This conjecture is based on the observation that 
the effect of skewedness (nondiagonality) should be negligible at large xi and small Qq from 
which the QCD evolution starts, i.e. one can neglect the effect of A at the input scale Qq. 
We refer the reader to Subsec. [IV A| where this question is investigated in detail by analyzing 



general tendencies of the QCD evolution in the x — Q-plane. 

Practical ways to estimate skewed parton densities through conventional parton densities 
were suggested in Ref. In particular, using symmetries of the kernels of the skewed 



QCD evolution, it was proposed in Ref. 0,^6| that the skewed parton densities at given x\ 
and x are expressed through the conventional parton densities computed at smaller x' = 
x\ — A/2. The numerical analysis of [47| showed that, indeed, it is more appropriate to 
compare the skewed parton densities to the usual parton densities taken at x' < x\. 

The factorization theorem for hard exclusive meson electroproduction off nucleon or 
nuclear targets enables one to experimentally measure the skewed parton distributions in 
the nucleon (nucleus) and the distribution of bare quarks within mesons and nucleons (the 
minimal Fock components in the wave functions of mesons and nucleons). The number 
of the skewed parton distributions is much larger than the number of parton distributions 
measured in inclusive DIS since the former includes, for example, transitions like N —>■ A, 
various distributions of spin, etc. Thus, it is reasonable to expect significant progress in 
the understanding of the structure of baryons and mesons from the investigation of hard 
diffractive processes. 
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III. HARD DIFFRACTIVE PROCESSES 



In this section, we discuss the implications of color coherence for hard diffractive processes 
on nucleons and nuclei. 



A. The phenomenon of Generalized Color Transparency 

The phenomenon of the complete disappearance of initial and final state interactions of 
spatially small-size wave packets of quarks and gluons, produced in hard exclusive processes, 
is usually referred to as the phenomenon of color transparency Initial theoretical ideas on 
the color transparency phenomenon were intuitive and based on model assumptions. For 
the history of the development of this subject and appropriate references, see, e.g., ]T). The 
recent progress in the theoretical understanding of the color transparency phenomenon is 
based on the systematic application of the QCD factorization theorem to hard exclusive and 
semi-exclusive reactions. In particular, the application of the factorization theorem to hard 
diffractive vector meson V production on nuclear and nucleon targets at sufficiently small x 
enables one to express the ratio of the forward differential cross sections for the nucleus and 
the nucleon as (compare to Eq. ([14])) 

da ~ / *+A-+V+A/ dt ( t = Q) 

da^+ N ^ v+N /dt{t = 0) ~ 

Eij Jo dz J dxxfj/Ajxi , x\ -x,t = 0,Q 2 ) H ij (Q 2 x 1 /x,Q 2 ,z)(pJ(z, Q 2 ) 
Eij Jo 1 dz J dx 1 f i/N {x x ,x 1 -x,t = 0,Q 2 ) H ij (Q 2 x 1 /x,Q 2 ,z) (j)J(z,Q 2 ) 



(15) 



where fi/A(xi,xi — %, t, Q 2 ) is the nuclear skewed parton density of parton i. 

The phenomenon, described by Eq. ([L^), can be called generalized color transparency 
PlJIlfnH . Although it resembles color transparency, the two phenomena lead to qualitatively 
different patterns of A dependence for sufficiently small x. Indeed, in case of the color 
transparency phenomenon, the complete disappearance of the initial and final state inter- 
actions of the qq pair, when x decreases, is predicted [J48|| . On the contrary, in the case of 



the generalized color transparency phenomenon, the initial state interactions with nuclear 
medium do not disappear but increase with a decrease of x. They are included into nuclear 
modifications of parton distributions (nuclear shadowing of parton distributions) the nuclear 
target f i / A (xi,x 1 - x,t,Q 2 ) (see, e.g., Fig. |). 

It is important to repeat that the considered example of diffractive production of vector 
mesons presents a particular case of the QCD factorization theorem which states that nuclear 
effects in hard exclusive processes at small x are calculable in QCD within the framework 
of the generalized (skewed) QCD evolution equation for the parton distributions of nuclei. 

Now let us consider two limiting cases of Eq. (13). In the limit of fixed x and Q 2 — > 



oo, the nuclear effects in the parton distributions disappear, i.e. fi/A( x i, x i ~ x ,tyQ 2 ) — > 
Afi/isr(xi,xi — x, t, Q 2 ). This leads to the more familiar phenomenon of color transparency 

da r + A^v + A /dt{t = 0) 

d<rr+ N -+ v + N /dt(t = 0) ' 1 ' 
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Eq. ([TBI) is a consequence of the disappearance of nuclear effects in the skewed nuclear 
parton densities as a result of the QCD evolution. The reason for such a disappearance is 
that the skewed (as well as conventional) parton densities at given x and Q 2 are determined 
by the input parton distributions at significantly larger x and smaller Qq where nuclear 
shadowing is unimportant. Thus, as one evolves from Q 2 , to Q 2 , the unshadowed large-x 
region contributes at smaller x and, consequently, diminishes the effect of shadowing. It 
explains how color transparency takes over generalized color transparency in the limit of 
fixed x and Q 2 — ■> oo. We also refer the reader to Subsec. |IV A| , where general trends of the 
QCD evolution in the x — Q-plane are discussed on the quantitative level. 

In order to illustrate the onset of color transparency numerically, we perform an analysis 
of the QCD evolution of the skewed and conventional nuclear parton densities. Details of the 
analysis are presented in Subsect. |1V B| . The limiting case, when x is fixed and Q 2 — *> oo, is 
presented in Figs. || and |7|. One can see that, at x = 10~ 3 and large Q 2 , the ratio of nuclear 
to nucleon skewed and conventional gluon parton densities increases and will eventually 
approach unity. This observation implies that the effect of nuclear shadowing vanishes and 
Eq. flTBD is justified in the discussed limit. 

In the other limit of fixed Q 2 and x — > 0, QCD predicts that nuclear effects begin to play 
an important role. This occurs due to the growing importance of nuclear shadowing at small 
x. This effect can also be interpreted as the decrease of the influence of color transparency. 
This limiting case for the conventional gluon parton densities is presented in Fig. ||. The 
ratio of the nuclear to nucleon conventional gluon parton densities becomes smaller at fixed 
Q 2 and x — > 0, which demonstrates the increase of nuclear shadowing. 

i.e. by the analysis of 



Within the models, inspired by the two-gluon exchange model p9 



the lowest-order diagrams of pQCD (see Ref. [48], for example), shadowing for the interaction 
of small dipoles is a higher twist effect which decreases oc 1/Q 2 when Q 2 increases. The 



difference between the predictions of |48[ and our analysis is due to the lack of Gribov 
diffusion to larger size configurations in the parton ladder in the lowest-order Feynman 
diagrams |IC|n|. Gribov diffusion has been observed long ago as a shrinkage of the diffractive 
peak in pp collisions with increasing energy. Another cause for the difference is that the 
QCD evolution equation predicts increasing parton densities at small x. This feature is also 
absent in the lowest-order Feynman diagrams. 



B. A brief review of experimental data 

The current subsection is concerned with hard exclusive reactions in the kinematic limit 
where color transparency plays a dominant role. Namely, these reactions are vector meson 
electroproduction off nuclear targets, coherent photoproduction of J /if), pion diffractive dijet 
production on nuclei and processes where the two-gluon exchange is prohibited. 

Exclusive meson production by longitudinal photons was discussed in the previous sub- 
section. Eq. ( fL5f) represents the ratio of the corresponding forward differential cross sections. 
In the limit of fixed small x and Q 2 —>■ oo, Eq. fllB]) has been proven in QCD within the 
leading a s log(Q 2 /A 2 ) logx approximation in |2lJ and within the leading a s log(Q 2 /A 2 ) 
approximations in M. 

Originally, in |7]j21[] , Eq. (|I~5|) was obtained using the conventional gluon parton densities. 
The factorization theorem for exclusive vector meson electroproduction requires that the 
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skewed parton densities are used instead of the conventional parton densities. However, 
regardless of whether skewed or conventional parton densities of the target are used, Eq. 
(|T5|) leads to Eq. ( |I6D because in the limit of fixed small x and Q 2 — > oo nuclear shadowing 
disappears in both skewed and conventional parton densities, as we explained in Subsec. 

An additional A dependence of the cross section for exclusive vector meson electro- 
production is also present. The cross section is proportional to the square of the nuclear 
form factor F\(t m i n ) = exp(—m 2 N x 2 3: jR 2 i /3), where —t min = m%x 2 B j is the minimal mo- 
mentum transfer to the target. However, due to the smallness of x B j considered here, 
XBj -C v^3/(-Ra^v) ~ 0.34/A 1 / 3 , the effect of the nuclear form factor is negligible in Eqs. 



(PD and flU). 

The behaviour, predicted by Eq. flT5|), has been verified experimentally in a number of 



phenomena. The E665 data on exclusive p meson electroproduction off nuclear targets f49 



reveals that the A dependence of the cross section increases with the increase of Q 2 . Such a 
behaviour is a characteristic feature of the color transparency phenomenon. However, this 
result should be considered with a grain of salt because part of the Q 2 dependence can be 
attributed to the fact that in the kinematics, covered by the E665 data, the coherence length 
l c is comparable with the nuclear radius and that within the data sample l c depends on Q 2 . 

Another example of processes, where color coherence plays a principal role, is coherent 
photoproduction of J/ip. The data on the cross section of coherent photoproduction of J/ip 
off nuclei |5(J at (E 7 ) = 120 GeV and x e ff ~ 0.04, where the shadowing effects in the 



gluon density are negligible, demonstrates that the cross section is proportional to A n with 
n m 1.4 ± 0.06 ± 0.04. This result is consistent with the idea in plj that J/ip, produced 
in exclusive photoproduction processes, originates from conversion of the photon to a cc 
quark pair at small relative distances, which hence interacts with a small interaction cross 
section. At the same time, the genuine J/^-nucleon cross section could be significantly 
larger because, in case of J/^-nucleon interaction, the cc quarks are at average distances 
within J /if) |T]|)2j. 

Hard diffractive electroproduction of vector mesons, described above, is dominated by 
small-size quark-gluon configurations, enforced by choosing the special kinematic condition 
of Q 2 — > oo. This condition selects small transverse size qq configurations in longitudinal 
photons. Another way to select small transverse size configurations, is to choose a special 
final state. An example of such a process is coherent diffractive production of mini-jets 
by pions in the reaction it + A — > jet(k t ) + jet(—k t ) + A. If we impose the restriction 
that only high transverse momentum k t jets are produced and that only the diffractive 
peak is investigated such that the nucleus A is intact, then it would insure only small qq 
configurations in the wave function of the incoming pions. 

The recent Fermilab experiment E791 data on diffractive dijet production by a pion 
projectile with momentum of 500 GeV/c off nuclear targets of carbon and platinum |53j is 
consistent with the predictions of the generalized color transparency |22| . The analysis of the 
E791 data demonstrates that cr(?r + A -> 2jets + A) oc A L61±a08 for 1.5 GeV/c < k t < 2.0 
GeV/c and a oc A im±am for 2.0 GeV/c < k t ||23|| . These results are consistent with the 
predictions based on color transparency. 

It is important to stress that this A dependence is significantly faster than the expected 
A dependence for usual soft diffraction discussed in Sect. M. Namely, the ratio of the cross 
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sections for platinum and carbon, ap t /ac observed in the E791 experiment, is approximately 
seven times larger than this ratio, estimated within the framework of soft diffraction, where 
a oc A 2 '\ 

In addition to the A dependence, the E791 experiment confirmed other expectations for 
the color transparency regime: i) The pion wave function at Q 2 ~ 10 GeV 2 (k t > 1.5 GeV/c) 
is described well by the asymptotic wave function calculable in pQCD, (^(z)! 2 oc z 2 (l — z) 2 , 
where z is the fraction of the longitudinal momentum of the pion carried by a quark in the 
infinite momentum frame [ 54 ,j5l)|| ; ii) The k t dependence of the cross section at k t > 1.5 



GeV/c is also consistent with the pQCD calculation of 

It is important to note that the observed A dependence of the cross section of diffractive 
dijet production is close to that found for exclusive photoproduction of J/ip |50| . This 
observation is in line with the QCD factorization theorem of [p] since the cross sections for 
diffractive dijet production and exclusive photoproduction of J /if) are expressed through the 
same nuclear skewed gluon distribution at similar light-cone gluon momenta. 

More experiments are needed to understand the dynamics of hard diffractive exclusive 
processes more clearly. We suggest measuring diffractive dijet production off nuclei by 
different projectiles, such as kaons, or to measure charmed jet production by real photons. 

Eqs. flTB] ) and ( |TED are very general. They have been proven using the QCD factorization 
theorem and the QCD evolution equation. According to the QCD factorization theorem for 
certain hard processes dominated by the scattering of compact system from hadron (nuclear) 
targets, the cross sections are expressed through the target parton (gluon) densities. But, 
as it follows from the QCD evolution equation, nuclear shadowing vanishes at fixed x with 
the increase of Q 2 . Hence, in this limit, Eq. (^) is justified. A similar prediction is also 
applicable to the processes where the exchange with vacuum quantum numbers in t channel 
is forbidden and where as a result the two- gluon exchange does not contribute. Below we 
will give a few examples of such processes. 

The QCD factorization theorem predicts the color transparency phenomenon for the zero- 
angle charge-exchange coherent processes, initiated by longitudinal photons in the limit of 
Q 2 — > oo and l/(2m^x) ^> Ra (so that the effect of t m i n in the nuclear form factor is small) 
but l e Jf <C r NN , where r N N is the average internucleon distance in the nucleus || 

da(Y L + A^p- + A')/dt = \(A\T + \A')\ 2 

da(7* L + B -»■ p~ + B')/dt \{B\T + \B')\ 2 ' 1 ' 

where A 1 (B') is the isobaric state of a nucleus with the atomic number A (B) with a neutron 
substituted by a proton; T + is the isospin-raising operator. 

In order to diminish the sensitivity of the discussed processes to the nuclear wave function 
of the isobaric states, it can be compared with the zero-angle charge-exchange coherent 
processes initiated by high-energy neutrinos or anti-neutrinos: 

dajvi + A^l + A^/dt = \{A\T + \A')\ 2 

da(p t + B ->l + B')/dt | (B\T+\B f ) | 2 ' 1 j 

where v\ is the neutrino corresponding to the lepton I. 

Another example of zero-angle charge-exchange coherent processes, where color trans- 
parency dominates, is the production of electrically charged pseudo-scalar mesons ||. In 



17 



the limit of Q 2 — > oo and l c > Ra but l e J^ <C r^Ar, QCD predicts for the reactions, initiated 
by virtual photons 

da(Y L + A^*- + A')/dt = \{A\T + \A')\ 2 (iq) 
da(Y L + B -> 7T- + £')M |(B|T+|5')| 2 ' 

It is important to emphasize that in the processes, described by Eqs. QT7f) , (|T8|), and (|19|), 
exchange by quarks is required. 

Similar predictions are applicable for the coherent production of ir + , p + and / mesons 
and pairs of mesons. 

The QCD factorization theorem || can be generalized to semi-inclusive processes off 
nuclei with large rapidity gap between the photon and target fragmentation regions. For 
zero-angle scattering, at Q 2 — > oo and fixed x, the onset of color transparency is expected 

^-ah*r +A^M + Rap. Gap + X') = A^-ah* r +N -> M + Rap. Gap + X') . (20) 
at at 

It is necessary to stress that in practical applications taking account of the leading twist 
effects in Eqs. (fl5|) , (|T6|) , and ( p0|) is insufficient since the small transverse size wave packet 
produced tends to expand to a normal hadronic size. In order to suppress this effect, one 
needs to impose the specific kinematic conditions - rather large Q 2 or sufficiently small x - 
so that the expansion will happen outside of the nucleus. 



C. Hard diffractive processes at HERA 

In this subsection, experimental results, obtained at HERA, are reviewed from the point 
of view of color coherence. 

Hard diffractive electroproduction of vector mesons on hydrogen 7* + p — > V + p, where 
V = p, 0, or J/ 1 ®, has recently been observed at HERA [jnj]. The data shows several 
phenomena, expected in QCD. Firstly, QCD predicts the increase of the cross section of 
hard diffractive processes, initiated by a projectile in a spatially small configuration of the 
size 1/Q 2 , as a power of the invariant energy of the collision^ 

cr(hard diffractive) w (xG p (x, Q 2 )^ . (21) 
Such a behaviour has been proven within the leading a s logx approximation in 0,0] and 



within the leading ot s lnQ 2 /A 2 approximation in |5|,|7|,^]. The result of Eq. (|2T| ) is in broad 
agreement with the energy dependence of cross sections for diffractive electroproduction 
of vector mesons p, <p, and J/ip observed at HERA \TA\. The rapid increase of the cross 
sections of hard diffractive processes with increasing energy is one of the signatures of the 
color transparency phenomenon |2~I]. 



4 Of course, one should use skewed parton densities in Eq. (21). However, in the kinematics of 
the experiments, the numerical difference between the skewed and conventional parton densities is 
rather small. 
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Secondly, the data have confirmed other QCD predictions for hard exclusive diffractive 
electroproduction of vector mesons such as the enhancement of heavy flavor vector mesons 
production (in contrast to moderate Q 2 processes, in which production of heavy flavors is 
suppressed), dominance of the longitudinal over transverse cross sections at large Q 2 , the 
universal dependence of cross sections of hard diffractive processes on the momentum transfer 
t, large cross sections of diffractive production of excited states of vector mesons, restoration 
of the SU(3) symmetry prediction for the ratio of cross sections for different vector mesons in 
the final state, etc. For the review and references see |I|flJT^|. The current data is consistent 
with the expectation of the onset of the hard QCD regime for hard diffractive processes at 



large Q Hence, the discovery of a variety of new calculable in QCD hard phenomena 



is expected. 

An important remark here is in order. All the above discussed hard exclusive diffractive 
processes were initiated by longitudinal photons. In such processes, the hardness of the 
photons (when the photons have large Q 2 ) leads to color transparency. This logic is inappli- 
cable to the processes, dominated by transversely polarized photons. In case of transverse 
photons, the soft component of the wave function (large size) may be as important as the 
hard component (small size). However, as a result of the Sudakov-type form factors, con- 
tributions of the soft component are suppressed at very large Q 2 . The relative importance 
of small and large dipole sizes in exclusive diffraction, initiated by transversely polarized 
photons, requires further studies. 

The importance of the soft component of the transversely polarized virtual photons 
manifests itself in the discovery of significant cross sections of diffraction in DIS |14| and in 



the leading twist nuclear shadowing in DIS [|1^,|I^]. For these processes, the hardness does 
not necessarily mean color transparency. The recent data of HERA [14] suggests that, for 
Q 2 > 6 GeV 2 , exclusive electroproduction of vector mesons is dominated by the compact 
configuration of the photon wave function. If so, generalized color transparency will be 
applicable to coherent diffractive electroproduction of transversely polarized vector mesons 
at these Q 2 . 



D. Hard diffractive electroproduction of photons 

The generalized vector dominance model (GVMD) is often considered as an alternative to 
the QCD explanation of the inclusive and diffractive DIS, for a recent review and references 
see [E7J. In the GVMD model, using the basis of hadronic states, the (virtual) photon 
interaction with the nucleon proceeds via the transition of the photon into different vector 
meson states, when the diagonal V + N — > V + N and non-diagonal V + N — > V + N 
transitions are accounted for. 

It was argued in [gg] that the GVMD models are successful in the theoretical descrip- 
tion of photoproduction and electroproduction at high energies since they possess (in an 
approximate way) basic features of theories with color coherence. 

However, it turns out that the process of diffractive electroproduction of photons in the 
DIS reaction 7*(Q 2 ) + p — > j*(QV) +p', which has been a subject of recent theoretical 
and experimental studies, can discriminate between the GVMD and pQCD pictures of DIS. 
Within the framework of pQCD, the ratio of the imaginary parts of the amplitudes of the 
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real and virtual photon electroproduction was found to be approximately independent of Q 2 
at small x [d^ 



i m A(Y(Q 2 )+P^i(Q 2 = Q)+p') „ 2 , 22) 
imA( 7 *(g 2 )+p^ 7 *(g 2 )+j9) ~ 1 ' 



The recent HERA data |Kj seems to confirm the result of Eq. (|22j). 

Within the framework of the GVMD model, which contains non-diagonal transition 
between the neighbouring vector mesons pl| , it was found in that 



ImA( 7 *(Q 2 ) + p -> 7 (Q 2 = 0) +p f ) 

imA(r(Q 2 )+p^r(Q 2 ))+p) 

(sr 6 M " V e \ /fsr e M 2 e Af£ x 

l i /„ M 2 + CP^fJ/ ^ /„ M 2 + Q^ nm f m Ml + Qi)- {Z6) 

The comparison of predictions of Eqs. (|22|) and (0) is presented in Fig. ||. One can readily 
see that the result of the GVMD model is consistent with the pQCD result only in a limited 
range of Q 2 , 1 < Q 2 < 5 4- 6 GeV 2 . Hence, the forthcoming HERA data on diffractive 
electroproduction of photons in the DIS reaction 7*(Q 2 ) +p — > 7*(Q 2 ) + p' in a wide range 
of Q 2 , would allow to establish the limits of applicability of the complementary description 
of DIS, which employs the GVMD models with elements of color transparency built-in. 



IV. NUMERICAL ANALYSIS OF THE DGLAP QCD EVOLUTION 

In this section, the onset of color transparency and generalized color transparency is 
studied numerically by performing the QCD evolution of parton densities of nuclei. Adopting 
a reasonable algorithm, we also analyze the paths in the x — Q-plane along which the QCD 
evolution proceeds. The knowledge of these paths enables one to substantiate the choice of 
the initial condition for the QCD evolution of the skewed parton densities and to explain 
the Q 2 evolution of nuclear parton densities. Also, the limits of applicability of the QCD 
evolution for nuclear parton densities are estimated. 



A. The QCD evolution paths 

The integro-differential DGLAP evolution equations, which govern the Q 2 and x depen- 
dence of parton densities, read as follows 




q s (y,Q 2 " 
g(y,Q 



2\ 



(24) 



In Eqs. fl24j) , q Nb (x,Q ) = qi(x, Q ) — qj(x,Q 1 ), with any i ^ j, is the quark flavor 
non-singlet distribution; q s (x,Q 2 ) = J2i{qi(x,Q 2 ) + (ji(x,Q 2 ) is the quark flavor singlet dis- 
tribution; g(x,Q 2 ) is the gluon distribution; a s (Q 2 ) is the running QCD coupling constant; 
P qq (x/y), P qg (x/y), P gq (x/y), and P gg (x/y) are the splitting functions calculable in pQCD. 
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Eqs. ( p4|) are often solved numerically. One example is the brute force method, when 
iterations in Q 2 are used. In particular, knowing the parton densities at Ql, one can find 
the parton densities at larger Q 2 using Eq. (|24]) as 



g» 3 {x, Q 2 ) = q NS (x, Ql) + Q2 n2 Q2 ° a 'f l) f % NS (y, Ql)P qq {-) • (25) 

vo Air Jx y y 

Similarly, Eq. (|24]) can be solved for the quark flavor singlet and gluon parton densities. By 
choosing the difference between Q 2 and Ql small enough and repeating the procedure of Eq. 
(p5p a desired number of times, each time using the parton densities from the previous step 
in Q 2 , one finds the solution of Eq. ([24]). From Eq. fl25f ) one can see that the parton densities 



at the input scale Ql determine the the parton densities at a larger scale Q 2 , Q 2 > Ql- 

The general trends of the QCD DGLAP Q 2 evolution are well known. As Q increases, 
the parton densities shift toward lower x due to the emission of soft partons. So, the QCD 
evolution proceeds along a path in the x — Q-plane, which extends from low Q and high 
x towards large Q and small x. The detailed knowledge of this path is very important. 
It enables one, for example, to estimate the influence of the input parton densities at the 
initial scale Qo on the result of the QCD evolution at higher Q and judge which regions 
of x at Qo contribute to the parton densities at Q > Q . Thus, our study of the general 
trends of the QCD evolution addresses the issues of the choice of the input distribution for 
the QCD evolution for the skewed parton densities and the Q 2 and x dependence of nuclear 
shadowing. It might also contribute to our understanding of the limits of applicability of 
the leading twist QCD evolution equation. 

In order to investigate a path in the x — Q-plane relevant for the QCD evolution, the 
following algorithm was adopted. 

• Fix the final point of the QCD evolution at some large Q and small x. We have studied 
three different final points at Q=5 GeV and x=10 -4 , 10~ 3 , and 10~ 2 . By solving Eq. 
(p5|) numerically, we find the parton densities at these final points. 



• Begin increasing the lower limit of integration in Eq. ( [25D from x to some x'. It will 
decrease the integration interval over y at the initial scale Qo, which contributes to 
the parton densities at the final points - the values of the parton densities at the final 
points will decrease. Choose x' such that the integration from x' to 1 in Eq. fl25|) 
gives the parton densities at the final points, which are half the parton densities, when 
the integration takes place from x to 1. Then, for each final point, i.e. for each pair 
(Q > Q ,x), x' represents the average x at Q . Take those x' as the initial point of 
the QCD path. 

• Now that the initial and final points of the path are known, one needs to connect them, 
i.e. one needs to find pairs {Qo < Q < 5 GeV, x") through which the path passes. At 
each Qo < Q < 5 GeV, x" is found by requiring that, when the integration in Eq. ( |25|) 
is from x' to 1, the parton densities at x" are half the parton densities at x", when the 
integration takes place from x to 1. 

The numerical solution of Eq. ([25]) was found using the CTEQ QCD evolution package. In 
order to illustrate the role, played by the shape of the input parton densities, two different 
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input parameterizations of Ref. [^2j were considered: CTEQ4LQ with the starting scale 



of the evolution Q = 0.7 GeV and CTEQ4D with Q =l.6 GeV. The parameterization 
CTEQ4LQ provides "valence-type" gluon and valence quark initial distributions localized 
around x = 0.2; the input sea quark distribution is quite small. CTEQ4D gives slightly 
rising input parton distributions as x decreases. 

The results for the path in the x — Q-plane for the gluon, u and s quark parton densities 
are presented in Fig. [TD|. CTEQ4LQ was used as an input. Each line represents a path, 
which begins at Qq—1.Q GeV and certain x' and ends at Q=5 GeV and x = 10~ 4 , or 10~ 3 , or 
10 -2 . Each curve passes through 1.6 < Q < 5 GeV and x", defined by the above algorithm. 
In Fig. [10], the upper panel corresponds to the implementation of our algorithm for the gluon 
parton density; the middle panel is for the u quarks; the lower one is for the s quark parton 
density. 

From Fig. |l0|one can see that the QCD evolution for gluons and quarks is different. The 
gluon density evolves in both Q and x; the u quark evolves mostly in Q and turns sharply 
to reach the final point only at the end; the trajectory for the s quark is a mixture of the 
trajectories for the gluons and u quarks. 

Note that the cross-over of the paths for the u quark in Fig. [10] is due to the following 
reason. At small x, x = 10 -4 4- 10 -3 , the u quark parton density is rising and large. Hence, 
its support could come only from x ~ 0.2, where the input parton densities are peaked. 
This explains why the initial points of the corresponding paths are located at large x. On 
the other hand, the u quark parton density has a local minimum at x ~ 10~ 2 and, thus, the 
large support region in Eq. (|25|) is not required. This leads to the initial point at x ~ 0.06 
for the path, which ends at Q=5 GeV and x = 10~ 2 . 

Fig. 10 presents a quantitative approach to the path of the QCD evolution. One could 
arrive at the same qualitative conclusions simply by examining the corresponding parton 
densities at different Q clS 8b function of x. 

In the sense of the proposed algorithm, Fig. |TD] answers the question: In Eq. (^), which 
regions of y at the input scale Qo are important at higher Q > Qq! For example, one can 
see from Fig. 10 that at Q=5 GeV and x = 10 -4 , or x = 10~ 3 , or x = 10 -2 half of the 



gluon parton density has come from x > 0.18 -4- 0.22 at Q =0.7 GeV. This observation can 
be applied to an analysis of the QCD evolution of skewed parton densities in the proton and 
conventional parton densities in nuclei. 

Let us assume that in the kinematics of a certain exclusive process, the skewedness 
parameter A = xsj is small. Then, its effect is negligible at large values of X\. We refer 
the reader to Subsect. |IID| for the definition of the kinematical variables. Since half of the 
contribution at small X\ and large Q has come from x\ > 0.2 at the input scale Qo, the 
influence of the effects of A will be decreased by at least 50% at small x\ and large Q. 
Therefore, this analysis justifies our claim that it is a good approximation to neglect the 
skewedness parameter A in the input for the skewed QCD evolution. 

Applying the result of Fig. |l^ to the QCD evolution of conventional nuclear parton 
densities, one can readily explain the decrease of nuclear shadowing with the increase of Q. 
Since the major contribution at large Q and small x comes from larger x at the input scale 
Qo where the effects of nuclear shadowing are not so significant, one would expect the effect 
of nuclear shadowing to decrease as Q increases. 

Fig. ITU] presents average paths. One can go one step further and study the influence of 
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various regions of y in Eq. (|25|) on the resulting parton densities in the final points and not 
only the average y that we called x' and x" . Such a study will help to answer the questions 
about the input for the skewed QCD evolution equation and the Q 2 behaviour of nuclear 
shadowing in more detail. 

To this end, let us introduce functions f g (x,z) for gluons, defined as 

£ i \ xG(x,Q = 5) 

f ^ z) = xG(x,Q = 5) z > (26) 

and similar functions for quarks and antiquarks. Here xG(x, Q = 5) is the gluon density at 
Q = 5 GeV and x = 10~ 4 , or 10~ 3 , or 1CT 2 ; xG(x, Q = 5) 2 is the gluon density at the same 
Q and x but originated from Eq. fl25|) , where the lower limit of integration is z. Thus, by 
definition f g (x, x)=l and f g {x, 1)=0. 

f(x, z) for gluon, u and s quark parton densities as a function of z are presented in Fig. 
|TT| . It helps to quantify further the importance of various regions of integration over y in 



Eq. (|24|). For example, from the upper panel of Fig. [11], one can see that at Q=5 GeV and 
x=10~ 4 , 80% of the gluon density has come from y > 0.1 at Q =0.7 GeV. 

We have also investigated how the QCD evolution path and the functions f g , u ,s( 

00 . Z^j £1X6 

affected by the choice of the input distribution for the QCD evolution equation. Figs. O 



and [13] represent the QCD evolution paths and functions f gi u,s{%, z) for the case, when the 
parameterization CTEQ4D was used as the input. One can see that the fact that CTEQ4D 
input parton densities increase slightly as x decreases, affects the QCD evolution paths. 
While the general trend remains unchanged, the initial points for the paths shift to lower x. 
In Fig. O, one can see that at Q=5 GeV and x=10~ 4 , 70% of the gluon density has 



come from y > 10~ 3 at Qq=1.6 GeV. Thus, if the skewedness parameter A is of order 10 -4 , 
its effect is negligible at X\ ~ 10~ 3 . Therefore, Fig. quantitatively demonstrates that A 
can be safely neglected in the input for the skewed QCD evolution. 

The comparison of Figs. |l^ and [l^ demonstrates that the proposed algorithm for finding 
the path, along which the QCD evolution for the gluon and s quark parton densities proceeds, 
does not depend on the initial scale Qq. Namely, the corresponding curves in Fig. ^ at Q=1.Q 
GeV pass through almost exactly the same x' as in Fig. ^ at Q =1.6 GeV. Thus, Figs. |TI] 



and [12] can be used to analyze the QCD evolution of the gluons and s quarks between any 
two scales. A radically different input for the valence quarks for the two parameterizations, 
considered here, makes the corresponding evolution paths in the region 1.6 < Q < 5 GeV 
very different, and, thus, the universality discussed above does not hold. 

In conclusion of this subsection, we would like to stress that our studies of the typical 
paths of the QCD evolution in the x — Q-plane presented in Figs. [T^ and [12] demonstrate that 
the input for the skewed QCD evolution can be taken as the input for the QCD evolution of 
conventional parton densities as long the input densities are constant or grow insignificantly 
as x —>■ at the initial evolution scale Qo- Such a behaviour at low x is constrained by the 
slow rise of cross sections for soft hadron-hadron interactions as energy increases. 

B. Numerical analysis of nuclear shadowing 

In this subsection we will present our numerical studies of the Q 2 and A dependence 
of skewed and ordinary parton densities in nuclei and nucleons using the QCD DGLAP 
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evolution equation. The studies serve the purpose of estimating numerically the regime of 
the onset of color transparency for hard exclusive and inclusive processes. 

As explained in Sec. |TJ, the kinematics of exclusive reactions dictates the use of skewed 
parton densities in evaluating cross sections of hard exclusive processes. In the limit of small 
and fixed x and large Q 2 , the onset of color transparency manifests itself as the absence of 
nuclear shadowing in the ratio 

o ( a n\ xg A (x 1 ,x 2 ,Q 2 ) , , 



As discussed in Subsect. |l 1B|, an additional A dependence in Eq. (£7]) due to the nuclear 
form factor is negligible at small Xsj ■ Hence, the contribution due to the nuclear form factor 
was omitted in our numerical analysis^. 

The numerical analysis of the conventional (non-skewed) QCD evolution equation was 
performed with the next-to-leading order in a s log(Q 2 /A 2 ) accuracy (NLO) using the 
CTEQ QCD evolution package. For the nucleon parton densities, we used the CTEQ4D 
parametrization with the starting evolution point Q$=1.Q GeV. The choice was motivated 
by the observation that the CTEQ4D gives a reasonable nonperturbative input: the gluon 
parton density is a rather flat function of Bjorken x at small x in accordance with physical 
intuition based on the slow rise with energy of soft hadron-hadron cross sections. 

In order to carry out the numerical analysis of the skewed QCD evolution equation, we 



used a modified version of the CTEQ QCD evolution package In Subsect. [TV A| , it 



is explained why the input for the skewed QCD evolution can be taken as being equal to 
the input for the QCD evolution for conventional parton densities. Note that the detailed 
investigation of this issue showed that the input distribution should be fairly flat for this 
conjecture to be valid. Consequently, we started our QCD evolution with the input skewed 
parton densities g(xi, x 2 , Qq) = xiG(xi, Qq) where G(x\, Qq) is the conventional parton 
density. The accuracy of our calculations is of the leading order in a s log(Q 2 /A 2 ) (LO) 
because the skewed evolution kernels of the evolution package are implemented with the 
same accuracy. Hence, for the input nucleon parton densities we took the leading order 
in a s log(Q 2 /A 2 ) CTEQ4L parameterization with the starting evolution point Qq=1.Q GeV 

In order to account for significant nuclear shadowing at 10 5 -j- 10 4 < x < 0.05 for gluon 
and sea quark nuclear parton densities and slight enhancement at 0.05 < x < 0.2 of the gluon 



nuclear parton density, we used the parameterization of Ref . [27 ] . This parameterization is 
based on the recent analysis of diffractive ep data of HERA, which revealed that the amount 
of shadowing in the gluon channel is approximately 3 times larger than was thought before. 

Since gluons and valence quarks evolve separately at small x, one does not need to 
model shadowing and antishadowing for the valence quarks when studying the gluon parton 
density. Should one need to study nuclear modifications of valence quarks, one can assume 
that the valence quark parton densities are shadowed at small x and then enhanced at larger 



Note also that for semi-inclusive processes when a sum over nuclear fragmentation processes is 
taken, there is no suppression due to the nuclear form factor. 
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x. However, one has to satisfy the baryon and momentum sum rules for the nuclear parton 
densities. 

Results of the QCD evolution for nuclei with A =12, 40, 100 and 200 are presented in 
Figs. ||-]8|. Figs. H and || are the results of the QCD evolution of the conventional (diagonal) 
parton densities. The conventional QCD evolution is relevant for hard inclusive processes 
off nuclei such as DIS. 

Fig. ^| illustrates decreasing nuclear shadowing due to the QCD evolution at high Q 2 and 
fixed x. The ratio of the nuclear to nucleon gluon parton density per nucleon Rg{%, Q) = 
Ga(x, Q)/AGn(x, Q) is plotted as a function of Q at x = 10~ 3 . One can see that the amount 
of shadowing for gluons, given by 1 — R G (x,Q), decreases by a factor of 5 as Q increases 
from 1.6 to 15 GeV for nuclei with A = 12 and A = 40 and by approximately a factor of 4.5 
for nuclei with A = 100 and A = 200. 

The onset of color transparency for exclusive processes (see also Eqs. (|i~6"D and (|27|)) 
is presented in Fig. [7]. The ratio of nuclear to nucleon skewed gluon parton densities 
Rg(xi, A, Q) = g^(xi,X2, Q)/Ag N (xi,X2, Q) is plotted as a function of Q at fixed X\ = 10~ 3 
for different values of x\ — X2 = XBj = A. One can see that the curves for different A for 
each A are remarkably close to each other. This is an interesting result because, even though 
the skewed gluon parton density is significantly larger than the diagonal one at small x, the 
effect of the skewedness A (asymmetry) cancels almost exactly in Rq{x\ 1 A, Q). It can be 
explained by the fact that the QCD evolution does not change the difference X\ — X2, which 
is fixed by the kinematics of exclusive reactions. From Fig. [7| one can see that the amount 
of nuclear shadowing for gluons, given by 1 — Rg(xi, A, Q), decreases by approximately a 
factor of 3.5 as Q increases from 1.6 to 15 GeV for nuclei with A = 12, A = 40, A = 100, 
and A = 200. 

Fig. |8| illustrates increasing nuclear shadowing in the limit x — > 0. At Q = 5 GeV, when 
nuclear shadowing is still significant, the ratio 1 — Rg{%, Q) increases by a factor of 10, when 
x is decreased from x = 10~ 2 to x — 10~ 4 for nuclei with A = 12, by approximately a factor 
of 8 for nuclei with A = 40, and by approximately a factor of 6 for nuclei with A = 100 and 
A = 200. 

Figs. ||, |8|, and [7| demonstrate quantitatively that, as Q increases, nuclear shadowing 
decreases, and the regime of the generalized color transparency turns into the regime of 
color transparency. 

C. The unitarity boundary for the inelastic cross section and the phenomenon of 

perturbative color opacity 

In this review we discuss coherent phenomena with nuclei at small Bjorken x. The main 
tool for the numerical analysis is the DGLAP QCD evolution equation. Then, the question 
arises: How low in x can one go before the evolution equation becomes unreliable? 

One way to approach the answer to this question was presented in Ref. @. It is based 
on the observation that the inelastic small-size qq dipole-nucleon cross section of Eq. (|1~3"|). 
being of the leading twist, dominates over the elastic gg-nucleon cross section, which is of 
the next-to-leading twist 

a inel > a el . (28) 
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When a mel ~ a el , the higher twist expression becomes of the same order as the leading twist 
term. Hence, the QCD evolution equation becomes inconsistent since it has generated large 
higher twist terms, which have been neglected in the derivation of the evolution equation. 

This argument demonstrates that a inel given by Eq. cannot grow infinitely at small 
x since that would destroy self-consistency of the DGLAP QCD evolution equation. The 
upper limit for a mel can be obtained from the constraint imposed by the unitarity of the 
S'-matrix for purely QCD interactions and the optical theorem M 
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(29) 



where r7v=0.8 fm is the r.m.s. radius of the nucleon, which is measured in the electropro- 
duction of vector mesons at large Q 2 via the t-slope of the differential cross section; (3 ~ 0.4 
is the ratio of the real to imaginary parts of the elastic scattering amplitude. 

In case of inclusive DIS on nuclear targets, the upper limit for the inelastic gg-nucleus 



cross section a A el follows from the unitarity restriction 
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where a A is the total gg-nucleus cross section. Eq. 



nuclei when a A 



2ttR a , where R A is the r.m.s. 



is applicable for sufficiently heavy 
radius of the nucleus. Using R A = l.lA 1 ^ 3 



fm, one can write Eq. (|3~0"D in the form 
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We studied Eqs. (^D) and ([T!]) numerically using the recently developed model for o mel 
The straightforward generalization of Eq. (|13|) for the color triplet (qq fluctuation) and 
color octet (qqg fluctuation) dipole scattering cross sections on the nuclear target reads 
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x'G A (x',X/b 2 )] a s (X/b 2 ) , 
x'G A {x',X/b 2 )} a s {X/b 2 ) . 



(32) 



In our numerical analysis, the nuclear gluon parton density x'G A (x', A/6 2 ) was parameterized 
as the product of the factor describing shadowing at 10 -5 < x' < 0.05 and antishadowing 
at 0.05 -T- 0.1 < x' < 0.2 according to Ref. | 27fl , and the gluon parton density of the proton. 



We used the leading order parton parameterization CTEQ4L in Eqs. ( p9|) and (|31|) since 
the expressions for the inelastic gg-nucleon (gg(yf-nucleon) and gg-nucleus (gg^-nucleus) cross 
sections have been derived in the leading a s (Q 2 ) log(Q 2 /A 2 ) approximation. 

The solutions of Eqs. (|29| ) and (|3TD are presented in Figs. |T4| - |T7| . For each given Q e ff — 
y/X/b, the liming value of x', which is denoted x/ im , at which Eqs. ( p9|) and ( plf ) are still 
satisfied, was found. Thus, the kinematics regions allowed by the unitarity restrictions lie 
to the right of the corresponding curves. 

In Fig. |ll, we plot xu m as a function of Q 2 for nuclei with A=12, 40, 100 and 200, 
and for the nucleon. One can see that, for nuclei with v4=200, 100 and 40, the kinematics 
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constraints are more stringent than for the nucleon. Thus, the breakdown of the DGLAP 
QCD evolution equation occurs sooner, i.e. at larger x for a given Q, for heavy nuclei than 
for the proton. 

Note that, for heavy enough nuclei, one can select scattering at central impact parameters 
by taking pairs of nuclei like 208 Pb and 40 Ca. Since the nuclear thickness at small impact 
parameters is about a factor of 1.5 larger than the average nuclear thickness, and since the 
average nuclear thickness is oc A 1/3 the plots for A=40, 100 and 200 roughly correspond to 
scattering at the central impact parameters on nuclei with A' « A/(1.5) 3 . 

In Fig. |T|, we plot xu m as a function of Q 2 for nuclei with A=200. We present three 
scenarios of gluon shadowing in nuclei. The solid curve corresponds to nuclear shadowing 
parameterized as in f27j (it is the same curve as in Fig. |i~4]) . The dashed curve corresponds 
to a smaller amount of shadowing. In this case, gluons are shadowed as the sea quarks in 



27J . The dotted curve is for the situation, when nuclear effects are neglected. In all three 



cases, at each given Q, xu m is larger than the corresponding xu m for the proton target. 

Figs. [16] and [17] represent the unitarity boundary for the inelastic qqg (color octet) fluc- 
tuation scattering cross section on nuclear and nucleon targets, which is given by combining 



Eqs. ( plf ) and (0). While the trend is quite similar to Figs. [14] and [y| the unitarity bound- 
ary is much more stringent - Xu m is greater by more than the order of magnitude at the 
corresponding Q e ff- 

The phenomenon of a tnel being large at small transverse size b and small x is called color 
opacity phenomenon. The regime of color opacity corresponds to a small coupling constant 
a s (Q 2 ) and large parton densities. Figs. [14] and |l| demonstrate that it is more advantageous 
to study the physics of color opacity with heavy nuclei than with the proton. The latter 
observation can be explained by the fact that the density of partons per nucleon, even when 
one takes into account significant nuclear shadowing, is larger for heavy nuclei than for the 
proton. 

The regime of color opacity is a new field of QCD. It becomes accessible at larger x for 
heavy nuclei than for the proton target. Hence, it should be feasible to study color opacity 
at LHC MM, HERA with nuclear beams MM, and eRHIC. 



V. SOFT DIFFRACTIVE DISSOCIATION ON NUCLEI 

The present section is concerned with applications of the principles of color coherence 
to the domain of soft physics - diffractive dissociation of hadrons on nuclei and nuclear 
shadowing in hadron-nuclear collisions. Color or cross section fluctuations constitute the 
basis of a phenomenological description of all the available data on diffractive dissociation 
and nuclear shadowing on light and heavy nuclei. 



A. Diffractive dissociation as a manifestation of the coherence length and color 

fluctuations 

As discussed in the Introduction, one of the principles of color coherence phenomena 
in QCD is that an energetic hadronic projectile consists of coherent quark-gluon config- 
urations of very different spatial sizes. These quark-gluon configurations are called color 
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fluctuations. They can also be called cross section fluctuations since the configurations with 
different transverse sizes interact with the target with different cross sections. Cross section 
fluctuations become significant, when the coherence length l c is large. The importance of 
the large l c in diffractive dissociation of hadrons was first discovered in |>6 . 



It is important to note that one has to distinguish two phenomenological descriptions of 
strong interactions. In the description, which is more relevant for hard processes, the color 
fluctuations are quark-gluon configurations of small sizes, which are almost fixed. The masses 
of such fluctuations are undefined. The other description, which is natural for diffractive 
dissociation, operates with physical particles of fixed masses. Then, it is more appropriate 
to consider the distribution over transverse sizes. Using the principles of color coherence, 
one can naturally relate the two descriptions of strong interactions. 

Let us consider diffractive dissociation of a hadron with mass m and large momentum 
Pi a b on a hadronic target T in the laboratory reference frame. The incoming hadron can 
fluctuate into a hadronic state with mass M* due to the uncertainty principle. The coherence 
length l c corresponds to the time, during which the incoming hadron stays in the state with 
mass M*, 

^^(v^^-V 7 ^^ (33) 

From Eq. (|33|), one can immediately see that the coherence length grows with the energy of 
the incoming hadron (which is proportional to pi a b)- 

If the coherence length is greater than the typical size of the target 2Rt, 

lc > 2Rt , (34) 



then all the fluctuations, whose masses satisfy Eqs. ( pB] ) and (0), interact coherently with 
the entire target. This fact has three consequences. Firstly, as energy grows, more diffractive 
states of ever-increasing mass can be produced. Secondly, the large l c corresponds to small 
momentum transfers to the target. It means that the target (nucleon or nucleus) has a large 
probability to remain in its ground state while the interaction takes place. Moreover, if 
the target is a nucleus, the suppression of the scattering amplitude due to the nuclear form 
factor is small. The first and second consequences manifest themselves in a large probability 
of coherent diffractive processes on nuclei. Thirdly, the coherence of the fluctuations means 
that they can be considered as frozen during the interaction with the target. This enables one 
to introduce a function -P(cr), which gives the probability to find a hadronic (quark-gluon) 
configuration, which interacts with the target with the cross section a, in the energetic 
projectile. -P(er) can be most naturally introduced within the framework of the basis of 
states, which consists of configurations of definite scattering cross sections with the target 



The vector of state of an energetic incident hadron \^f) can be represented as a coherent 
superposition of eigenstates j^) of the scattering matrix 

|¥) = $>*!**> » (35) 

k 

where 
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ImT\V k ) = t k \* k ) , 
El c fc! 2 = 1 - ( 36 ) 

k 

Here, T is the scattering operator, and t k is the imaginary part of the eigenvalue corre- 
sponding to the eigenstate \^k)- 

Since various states \fy k ) interact with the target with different cross sections a k , which, 
by the optical theorem, are related to the imaginary part of the scattering amplitude t k , 

Ok = tk , (37) 

the coherent superposition of the eigenstates, which forms the final state emerging after the 
scattering, could be different from the initial state. Thus, the formalism of eigenstates of 
the scattering matrix is natural for describing diffractive dissociation of hadrons. 

Using Eqs. (|3~5)-(|3"7), diffractive dissociation can be presented as follows. Diffractive 



scattering occurs when the final state carries the same quantum numbers as the initial state, 
i.e. whenever the initial state overlaps any \^k)- Then, the total diffractive differential cross 
section at t = can be presented as 

rdp&ff = J_y mT |^|2 = J_yi \2 f 2 J_, 2x (3g) 

In Eq. (|3^), we have used the completeness of the set of states \^k) and the optical 
theorem (j37|). We have also introduced the second moment (a 2 ) of the distribution over 
cross sections -P(cr). Similarly, the elastic diffractive differential cross section at t = reads 

sdayi = J_ mi mT|*>| 2 = — (V| c k \ 2 t kY = —( a) 2 , (39) 

where (a) is the first moment of -P(cx). 

Subtracting the elastic cross section from the total diffractive cross section, one obtains 
the diffractive dissociation cross section (inelastic diffractive cross section) 

,dayei = sdoyiff _ .day = J_, _ x 



Eq. (f40| ) was first derived in |68| in order to describe diffractive dissociation within the 
framework of cross section fluctuations. Eq. ([|y) explicitly demonstrates that diffractive 
dissociation occurs only if different components of the incident hadron interact with 
the target with different strengths t k , i.e. when cross section fluctuations take place in the 
wave function of the incoming hadron. 



B. Properties of -P(cr) 

The formalism of scattering matrix eigenstates provides an economic way to take into 
account the composite structure of energetic projectiles. Due to the frozenness condition, 
each state scatters independently on the target and then contributes coherently to 
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the final state with the probability \ck\ 2 , see also Eq. fl35|). Then, one can introduce the 
distribution over cross sections P(&) Jl ] 

P(<r) = 'Z\c k \ 2 5((T-(T k ) . (41) 
k 

One can see from Eq. fl4"T| ) that it is the composite nature of hadrons that gives rise to the 
term cross section fluctuations. 

The cross section distribution function P(cr) comprises perturbative and non- 
perturbative information on the wave function of the hadron. The modelling of P(cr) involves 
pQCD calculations, valid at small a, and constraints on first moments of a, extracted from 
experiments. In addition, it is required that P{cr) — > at o — > oo. 



The quark counting rule fixes the low-cr behaviour of Ph(&) for hadrons [6^ 

P h (a) cx a n "~ 2 , (42) 

where n q is the number of valence quarks in the hadron. Thus, for protons, when a is much 
smaller than the average value of a, i.e. when o <C (a), one obtains that 

P p (a) ~ a . (43) 



For pions, the quark counting rule predicts at a <C (a) [|69 | 



P 7r ((r) ~ const . (44) 

However, one can go further than the simple counting rule in case of pions. In order to 
find P n (cr) at a <C (cr), one can apply the QCD factorization theorem and find the 

constant entering Eq. 



6/2 

P ^ a) = 5a s (4kf)xG N (x,4k 2 ) ' (45) 

In this equation, a s (Ak 2 ) is the QCD running coupling constant; Gat(x, 4fc 2 ) is the gluon 
distribution in the nucleon; x = 4k 2 /s n N , where s^n is the center of mass energy squared; 
k\ (x l/b 2 , where b is the transverse size of the qq pair of the wave function of the pion; f n 
is the constant for the 7r — > {iv decay. 

Note that the calculations of the cross section of exclusive electroproduction of vector 
f2~l"H and that of P n (cr) use a similar technique. Therefore, the confirmation of 



mesons 



the results of by the ZEUS data at HERA [|1J] provides an indirect confirmation of the 
calculation of P 7r (a (a)). 

Using the QCD factorization theorem, one can compute -P 7 (cr) for photons at a (a). 
It was found that the behaviour of -P 7 (cr) at small a is divergent 



PJa) oc - . (46) 



Such a behaviour is consistent with the quark counting rule of Eq. (^Sp. Since the quark 
and antiquark in the photon wave function are created at the same point, they have a large 
probability to stay close to each other. 
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Suggested parameterizations of the proton and pion distribution functions P(o) describe 
correctly the above discussed behaviour at small and large a [B{| 



PJa) = N(a,n) "J^ e -(^°)"/(£M>)" 
a/cr + a 

P^a) = N(a,n)e- {r7 - aoyi/{nao)n . 



(47) 



Pp(a) and P n (a) as functions of a are presented in Fig. [18] |p9| . 

The free parameters of the parameterizations (fi7|) are chosen in order to satisfy the 
following constraints on the first few moments of the distribution P(cr): 



JP(a)da=l , 
!P{a)ada=(a) , 
J P(a)a 2 da = (a 2 ) = 
I P(a)(a - (a)) 3 da 



(1 
i 



(48) 



The first equation is the normalization of P{cr). The first moment (a) is the total hadron- 
nucleon cross section. 

The third of Eqs. ( |4"8"| ) is a relation between the second moment (a 2 ) and the quantity 
uj a . Using Eq. (|40|) , one can write that 



, inel 



, dt) t=0 



el 



'da\ 
v dt ) t=0 



{° 2 ) - (qY 

{a) 2 



(49) 



which demonstrates that uo a is proportional to the amount of diffractive dissociation. uj a 
can be extracted from the experiment. From the data on the total cross section of neutron- 
deuteron scattering |70| and from the theoretical analysis of |f7I |, which used the proton- 
deuteron data ff2| , u; o -=0.25 for pi ab ^ 



200 GeV/c. Another source of information on uo a is 
the data on diffractive dissociation on hydrogen. The data leads to a similar value of io a 



The last equation in Eq. ([|8|) is based on the analysis of \W\ of the data on coherent 
diffractive dissociation of protons on deuteron [73[]. It indicates that P p (cr) is symmetrical 
around its average value. In order to impose further restriction of higher moments of P{&), 
one has to consider diffractive dissociation on heavier nuclei, for example, on 4 He. 



C. Total and diffractive dissociation cross sections of hadron-deuteron scattering 

In the previous subsection we demonstrated that in order to study the second moment 
of -P(cr), (a 2 ), one can study either diffractive dissociation on hydrogen (see Eq. (^0|) ) or the 
total hadron-deuteron cross section. In this subsection, the latter approach is considered in 
detail. 



As discussed in Subsec. [II A| , there is a deep connection between high energy diffraction 
of hadrons on the proton and nuclear shadowing in the total hadron-nucleus cross section. 
Using the space-time picture of high energy scattering, which follows from the quantum field 
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theory, V. Gribov has derived the total cross section of hadron-deuteron scattering [25 
(see also Eq. @) 

id 2 a h+p^X+p da h+p^h+p^ 



hD hp , hn 

a tot — a tot ~r a tot 



where a t £ t and <t^ ™ are the total hadron-proton and hadron-neutron cross sections; 
d 2 a/dM 2 dt is the cross section for producing a diffractive state X with mass M in the 
reaction h + N — > X + N; t is the squared four-momentum transfer; S(t) is the electro- 
magnetic form factor of the deuteron. For simplicity, spin effects in the wave function of the 
deuteron and the real part of the diffractive scattering amplitude have been neglected. 

In spite of different space-time evolution pictures of scattering processes in quantum field 
theory and in non-relativistic quantum mechanics, Eq. ([50]) generalizes the earlier result of 
|74j , when only the second term in the bracket of Eq. ([5(]) (elastic shadowing) was taken 



1 r , n . ./ f , 9 d z a n+p ^ +p do n+p ^ n+p s , , 

* S{U ){ i li'- mqr + —; r -), (50) 



into account. 

Eq. (|50D relates the amount of nuclear shadowing in the total cross section of hadron- 
deuteron scattering to the cross section of diffractive dissociation of the hadron on the proton. 
Nuclear shadowing due to diffractive dissociation (inelastic diffraction) is called inelastic 
shadowing because it is caused by the excitation and propagation of inelastic intermediate 
states. Fig. [1] depicts that the inelastic shadowing term appears due to the interference 
between the diffractive scattering of the incident hadron off the proton and the neutron 
of the deuterium target. Although inelastic shadowing corrections are rather small for 
hadron-nucleus collisions, inelastic shadowing is the only mechanism, which relates nuclear 
shadowing in inclusive DIS on nuclei to diffraction in DIS on proton. 

Within the formalism of cross section fluctuations, Eq. (|50D was rederived in |75|] , 



o% = 2(a)-!&Je-*f'S{q) t Pq. (51) 

Here, it is assumed that a v tot = cr" ot = (a) and that the spin effects in the deuteron wave 
function and the real part of the diffractive scattering amplitude are negligible. Eq. (|51~D 
explicitly demonstrates how (a 2 ) and uj a can be extracted from the shadowing term in the 
total hadron-deuteron cross section. Using the data on elastic and inelastic shadowing in 
a tat EOI an d viat [EH ( see a l so EH)> one finds that u a (p,n) = 0.25 for the average values of 



Piab, Kafe=200 GeV/c. 

For a pion projectile, the value of uv(7r)=0.4 at £^=200 GeV/c was extracted in |76 
from the data on ttD scattering of J72J. This value of u a is in accord with the theoreti- 
cal expectations, ensuing from the quark counting rule and the triple Pomeron exchange 
model. The larger value of uj a {^) for pions than for protons means that the distribution over 
cross sections P{cr) is broader for pions than that for protons. This fact is reflected in the 
parameterizations (|4"7|). 

A direct comparison of Eqs. (f4~9"D, (|50|) and ( |5T]) enables one to identify uj a as the ratio 
of the diffractive dissociation to elastic cross sections 

(a 2 ) -(a) 2 ( f f „ tl ,d 2 a h+p ^ x+p - , x /, r, ,da h+p ^ h+p 



{ay 



Eq. (p2|) explicitly demonstrates the connection between diffractive dissociation, inelastic 
nuclear shadowing, and cross section fluctuations. 
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The third moment (a 3 ) of Phip) can be found by considering diffractive dissociation of 
the hadron h on deuteron. The reaction p + D — > X + D was analyzed in [|69|, where the 
ratio of the proton coherent diffractive dissociation cross sections on deuterons and hydrogen 
was presented as 



/ d<Jdiff \P D // d(Jdiff \PP 
\ dt >t=J \ dt Jt=o 



1 (a 3 ) -(a) (a 2 ) 
2n 2 (a 2 ) - (a) 2 



d 2 p±S(p±) . 



(53) 



Using Eq. (|53f) and the data on coherent diffraction of protons on deuterium |73 
at J P p (cx)(<7 — (a)) 3 da ~ (see p9| for details). 



one arrives 



D. Diffractive dissociation on He 4 



In order to study further the importance of color fluctuations in energetic hadrons and 
higher moments of the distribution P(a), one should consider heavier nuclei. The only 
existing experimental data on inclusive coherent diffraction on light nuclei is the data on 
proton diffractive dissociation on 4 He in the reaction p + 4 He — > X + 4 He W7\ . 



Using the cross section (color) fluctuation formalism along with the Gribov-Glauber 
multiple scattering formalism p5| , [74]| , the ratio of diffractive differential cross section on 4 He 
and on the proton at t — can be presented as a series over moments of P p (cr) [78| 



pHe 




t=Q 



57 3 /(a 5 ) 



(a 2 )\ 59 7 2 /(a 4 ) 27 (a 2 ) 2 32 (a 3 ) 



(a) 2 J + 48^V(^) 4 
4 (a 3 ) (a 2 ) l(a 4 )\ 
5 



59 (aY 



5 (ay 



59 (a) 3 J 
(54) 



Here, (a) =41 mb is the total hadron-nucleon cross section at pi a b = 200 GeV/c. The 
parameter 7, defined as 7 = (a) / n(a + (3), emerges from the integration over transverse mo- 
menta of the proton-nucleon scattering amplitudes with the 4 He ground state wave function; 
ct=23 (GeV/c)~ 2 is the slope of the ground state wave function of 4 He, which is parame- 
terized in the exponential form /?=10.5 (GeV/c)~ 2 is the slope of hadron-nucleon cross 
section, which is taken as the average of the slopes of the elastic and diffractive cross sections 
|7§. This yields 7= 1.01. 



Then, using the parameterization (17) in order to compute the higher moments of P p (<j) 
Eq. (H) leads to " 



( da djff \ pHe ,/ d(Jdiff \ PP 

v dt )t=o'\ dt h=c 



6.8 -T- 7.6 . 



(55) 



The lower value of r corresponds to n = 2 in Eq. ( f47|) , and the upper value of r corresponds 
to n = 6 and n = 10 in Eq. (|47|) . 

Eq. (^) is in an excellent agreement with the experimental data on proton coherent 
diffractive dissociation on helium [77]. Integrating the data over the available range of 



diffractive masses 2.5 < M 2 < 8 GeV 2 , it was found in 1731 that 



(56) 
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Here, the main error originates from the extrapolation of the data to t = 0. It is important 
to note that, although the comparison to the theoretical value of r implies the integration 
over all diffractive masses M 2 , the correction, related to the inclusion of non- measured small 



diffractive masses, is small 78 



The remarkable agreement between the theory and experiment demonstrates that, in- 
deed, significant color fluctuations (cross section fluctuations) are present in energetic pro- 
tons. The phenomenological description of those fluctuations by Eq. fl4"7|), which is based 
on the notion of color fluctuations, successfully reproduces the existing experimental data 
on diffraction on helium and heavier nuclei (see Subsect. |V E) with an a accuracy of 10% 



without introducing any new free parameters. The latter fact is a sign of universality of the 
description that uses the distribution over cross sections -P(cx). 

One can also compare the theoretical prediction with the data on the total diffractive 
cross section of proton- 4 He scattering. Combining the Gribov-Glauber multiple scattering 
formalism, cross section fluctuations and the parameterizations fl4"7|), the total coherent 
diffractive cross section of proton- 4 He scattering o d ^J_^ was obtained in |79| 



°He-4 = 3-5 mb, n = 2 , 

a d ^ f J_ 4 = 3.7 mb, n = 6, 10 . (57) 
The result of Eq. (|57|) is in a good agreement with the experimental data of [77], inte- 



grated over 0.04 GeV 2 < \t\ < 0.4 GeV 2 and over diffractive masses GeV 2 < M 2 < 10 
GeV 2 : 

a% f /_ 4 (exp) = 2.6 ± 0.9 mb . (58) 

In Eq. fl58|), the error is dominated by the uncertainty of the extrapolation of the experimen- 
tal data into the unmeasured regions of small t < 0.04 GeV 2 and small diffractive masses 
M 2 < 1.2 GeV 2 . 

Fig. ITyi presents a comparison of the theoretical predictions for the total cross sections 



of coherent diffraction of protons and neutrons on nuclei to the available experimental data. 
The value of a^L^exp) , given by Eq. (|58D , is presented as a full circle. The theoretical 
prediction for cr^{^ 4 , given by Eq. fl5"7|) , is presented by two dashed curves that correspond 
to n = 2 (lower curve) and n = 6 and n = 10 (upper curve). 



E. Diffractive dissociation on heavy nuclei 

Cross section (color) fluctuation formalism can also be successfully applied to diffractive 
dissociation on heavy nuclei. For heavy nuclei, one can neglect the slope of the elementary 
hadron-nucleon elastic scattering amplitude as compared to the slope of the nuclear wave 
function as well as correlations in the nuclear wave function. Then, in this approximation, 
the coherent diffractive dissociation cross section o\ht of hadrons \h) on the nucleus with 



the atomic number A can be presented as [80 



hA 

a diff 




daP h (a) J2 I (h\F(a, b)\n)\ 2 -(J daP h (a) | (h\F(a, b) \h) \f) , (59) 
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where b is the impact parameter; |n) is the diffractively produced final state; F(a, b) is the 
Glauber profile function, defined as F(a,b) = 1 — e~ aT ( b ^ 2 , where T(b) = /f^ paQ), z)dz, 
PA{b,z) is the nuclear density. Using completeness of the states \n), o\tff can ^ e presented 
as 

°*ff = J d% { J daP h (a)\(h\F 2 (a,b)\h)\ - ( / daP(a)\(h\F(a,b)\h)\) 2 ) . (60) 

The theoretical prediction of Eq. (|60"D for protons (neutrons) is compared to the experi- 
mental data in Fig. [TJ| where <?dtff * s presented as a function of A. The three solid curves 
correspond to the calculation with n—2, 6 and 10 of the parameterization (f47j). The A 
dependence of (Jdtff is °f the approximate form A ' 8 for A rs 16 and A 0A for A rs 200. There 
are only two sets of experimental data on coherent nuclear diffraction of nucleons on heavy 
nuclei - the data on n + A — > pn~ + A for the mass interval 1.35 < M < 1.45 GeV and 



the data for emulsion targets p2[. In Fig. O, the data of Hfl is presented as stars and the 



data of HH is presented as triangles. 



As one can see from Fig. [19], the theoretical prediction of the A dependence of coherent 
diffraction of protons and neutrons on nuclei, based on color fluctuations accumulated in 
P p (cr), reproduces the experimental data well. Note that the data points for the reaction 
n + A —>■ pn~ + A lie systematically below the theoretical (solid) curves because the data 
points represent only one channel (a particular final state) of the total diffractive cross 
section. 

For pion projectiles, the approach, which uses P n (a), also describes the available exper- 
imental data well [p0[] . 

Inelastic diffraction (diffractive dissociation) on nuclei arises due to scattering at impact 
parameters that are rather close to the nuclear surface, where the nucleus is not so opaque 
and the nuclear density is small. The dominant contribution to inelastic diffraction comes 
from the fluctuations of the projectile, whose size is close to the average size of the projectile. 
This fact can be seen from the approximate diffractive cross section <J^fj^ hA , which was 
computed by expanding F(a, b) about a = (a) in Eq. (E^) [K0| 

,., /^\2 

appr hA 



a diff 



Jd%T(b) 2 e~^ . (61) 

Eq. (|6l|) is quantitatively accurate for A < 50 and qualitatively good for all A. It demon- 
strates that the A dependence is mainly determined by the value of (a). For A > 50, Eq. 
(|Bl~|) gives a weaker A dependence than the exact formula (|59"1). The deviation at large A 
should come from configurations with relatively small a. 

It is important to note that the numerical analysis of |8(J shows that it is impossible to 
study the dominance of small size contributions in reactions of coherent nuclear diffraction 
simply by increasing A, because it requires unreasonably heavy nuclei. Consequently, one 
should explore the domain of hard diffraction in the search for the dominance of small-size 
configurations. 

F. Cross section fluctuations in heavy-ion collisions 

Color fluctuations also play an important role in inelastic reactions with nuclei. The effect 
of color fluctuations in the wave function of hadrons on multiplicity of binary nucleon-nucleon 
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collisions and the transverse energy production in central high-energy nuclear collisions was 
calculated in |||]. It was found that the color fluctuations account for the large multiplicity 
of binary nucleon-nucleon collisions and the transverse energy fluctuations found experimen- 
tally. The fluctuation parameter io a (see Eq. (f49[) ) increases at ^fs < 60 GeV and decreases 
for ■sfs > 400 GeV. Thus, u) a is likely to be close to its maximum, u a ~ 0.3 -v- 0.35, for the 
energies to be studied in AA collisions at RHIC. Therefore, the effects of color fluctuations 
in the hadron wave function are relevant for the RHIC program and should be included in 
the event generators of central nucleus-nucleus collisions. 

VI. CONCLUDING REMARKS 

We find that QCD unambiguously predicts the onset of color transparency in hard coher- 
ent meson electroproduction on nuclei at small Bjorken x and Q — > oo. Thus, coherent hard 
phenomena are now becoming an effective instrument in the investigation of the quark-gluon 
structure of hadrons and nuclei. 

At the moment, there are many uncertainties in the subject of extremely small Bjorken 
x phenomena, where the QCD evolution equation becomes inapplicable. However, the con- 
clusion that the onset of the new QCD regime of large parton densities can be achieved in 
eA collisions at HERA and pp collisions at LHC, seems to be model-independent. 
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FIGURES 




FIG. 1. Diagrams demonstrating the relationship between DIS diffraction on protons in the 
reaction 7* + p — ► X + p' and the leading contribution to nuclear shadowing in inclusive DIS on 
nuclei. 



40 




FIG. 2. The ratios of the nuclear and nucleon gluon parton densities Ga/AGn and the quark 
parton densities qA/Aqx for lead (Pb) and carbon (C) as a function of x. The dotted, dashed and 
solid curves correspond to Q=2, 5 and 10 GeV, respectively. 
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FIG. 3. l c rriNX as a function of Q at x=10 -4 , 10~ 3 , and 10~ 2 as given by Eqs. (|l0| ) and (|IT|), 
The solid curves correspond to the parameterization of the diffractive parton densities of [31] . The 
results, obtained with the parameterization, which also includes the low-/? tail (see explanation in 
the text), are given by the dashed curves. 
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FIG. 4. The block structure of the amplitude for hard exclusive production of mesons M by 
longitudinally polarized photons 7^ in the reaction 7£ + p — > M + B, where p is proton and B is 
a baryon. 
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FIG. 5. Diagrams demonstrating the relationship between the gluon induced hard diffraction 
on protons and the leading twist contribution to nuclear shadowing in exclusive hard vector meson 
production in DIS on nuclei. 
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FIG. 6. The ratio of the conventional nuclear to nucleon gluon parton densities Rg(x, Q) as a 
function of Q at x=lCT 3 . The input distribution for the QCD evolution is CTEQ4D. 
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FIG. 7. The ratio of the skewed nuclear to nucleoli gluon densities Rc(xi,X2,Q) as a function 
of Q at xi=l(T 3 and at A = 0.5 x 1(T 3 , A = 0.75 x 10" 3 and A = 0.9 x 10~ 3 . The three curves 
for each A are barely distinguishable. The input distribution for the QCD evolution is CTEQ4L. 
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FIG. 8. The ratio of the conventional nuclear to nucleon gluon parton densities Rg(x, Q) as a 
function of x at Q=5 GeV. The input distribution for the QCD evolution is CTEQ4D. 
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FIG. 9. The ratio of the imaginary parts of the amplitudes for real and virtual photon produc- 
tion by virtual photons as a function of Q 2 . The straight line is the pQCD prediction [59]. The 
result of the calculation using the generalized vector meson dominance model (see Eq. (^) is given 
by curve. 
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FIG. 10. The QCD evolution path in the x — Q-plane for the gluon, u quark, and s quark 
parton densities for the three final points, x = 1CT 4 , 10~ 3 , and 1CT 2 . CTEQ4LQ with Qo=0.7 
GeV is used as an input for the QCD evolution. 
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FIG. 11. Functions f(x,z) as a function of z and at x = 1CT 4 , 10 -3 , and for the gluon, 
u quark, and s quark parton densities. CTEQ4LQ with Qq=Q.7 GeV is used as an input for the 
QCD evolution. 
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FIG. 12. The QCD evolution path in the x — Q-plane for the gluon, u quark, and s quark 
parton densities for the three final points, x = 10~ 4 , 10~ 3 , and 1CT 2 . CTEQ4D with Q =1.6 GeV 
is used as an input for the QCD evolution. 
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FIG. 13. Functions f(x, z) as a function of z and at x = 10~ 4 , 10~ 3 , and 10~ 2 for the gluon, u 
quark, and s quark parton densities. CTEQ4D with Qq=1.6 GeV is used as an input for the QCD 
evolution. 
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FIG. 14. The unitarity boundary for the inelastic (/(/-nucleus cross sections for nuclei with A=12, 
40, 100, and 200. The unitarity boundary for the inelastic gg-nucleon cross section is presented as 
a thick solid curve. 
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FIG. 15. The unitarity boundary for the inelastic qq- nucleus cross sections for nuclei with 
^4=200. The amount of nuclear shadowing for gluons is varied: the solid curve corresponds to the 
full amount of shadowing, the dashed curve corresponds to the gluons shadowed as the sea quarks, 
and the dotted curve is for gluons without shadowing. 
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FIG. 16. The unitarity boundary for the inelastic qqg (color octet)-nucleus cross sections for 
nuclei with A=12, 40, 100, and 200. The unitarity boundary for the inelastic ggg-nucleon cross 
section is presented as a thick solid curve. 
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FIG. 17. The unitarity boundary for the inelastic qqg (color octet)-nucleus cross sections for 
nuclei with ^4=200. The amount of nuclear shadowing for gluons is varied: the solid curve corre- 
sponds to the full amount of shadowing, the dashed curve corresponds to the gluons shadowed as 
the sea quarks, and the dotted curve is for gluons without shadowing. 
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FIG. 18. The distribution over cross sections P(cr) for protons and pions, presented in Eq. 
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FIG. 19. The cross section of coherent diffraction dissociation of protons and neutrons on nuclei 
as a function of the atomic number A. The solid curves are the theoretical prediction of Eq. (|59|). 
The data on the reaction n + A — > pir~ + A [81] is presented as stars and the data for the emulsion 
targets [82] is presented as triangles. The data point, presented by the full circle, corresponds to 
A = 4 ( 4 He) and is extracted from [77] (see Eq. (|58|)). The theoretical prediction for coherent 
diffraction on 4 He is given by the dashed curves (see Eq. (|57|)). 
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